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ABSTRACT 


A three-region static 0(F) curve is calculated from two empirical 
hyperbolic functions that use five material parameters; very close agree- 
ment with experimental data is obtained in a practical range of applica- 
tion. After the limitations of a leg model are pointed out, the 
inelastic-switching parameter p of a tapered leg is calculated as a 
function of material properties, geometry, and flux level; this param- 
eter is useful for practical multipath cores whose legs usually vary in 
width. Based on a parabolic model for p(<fi) , average p of a constant- 
width leg and a tapered leg are calculated and plotted vs. flux change; 
these plots are applicable to circuits involving partial or full switch- 
ing, Flux- di vi s i on ratio D between two legs in parallel vs. MMF drive 
NI , applied to a third leg in series, is calculated in terms of leg 
parameters for a general case, and then simplified to special cases, 
e.g., unloaded core, or NI co . For relatively low N I , leg parameters 
are improved by considering the nonlinearity of the function peak 0 vs, 

F. For relatively high NI , the effect of leg dimensions on D is ex- 
pressed more directly. In view of the complexity of the analysis, the 
agreement with experimental data, especially in the no-load case, is 
quite good. Experiments were also performed in order to study flux 
division in a nonsaturable core. Modification of the parabolic switch- 
ing model to account for variation of the dynamic threshold with flux 
for low MMF results in closer agreement with experimental 0(0) data. 
Problems in determination of leg dimensions are discussed relative to 
two cores that were designed for f 1 ux- di vi s i on experiments. A new model 
is proposed for the process of slow flux switching in a core having a 
re-entrant shape. Laboratory techniques for determination of switching 
parameters are presented, and examples are given. The unsetting effect 
is measured as a function of applied MMF; the results are in full agree- 
ment with the qualitative analysis in a previous report. 
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PREFACE 


This is the second report on flux switching in multipath cores. In 
the first report, all available published and unpublished material re- 
lated to this subject was organized and evaluated, and new material was 
added. Conclusions were drawn about areas that needed further investi- 
gation. An attempt is made here to fill some of these gaps. Since this 
report is a direct continuation of the first report, reference to Report 1 
is made frequently in order to save space. The reader is advised to 
have Report 1 on hand; otherwise the chain of reasoning in many occasions 
will be hard (or even impossible) to follow. 

The material in this report is divided into two sections: The first 

section contains theoretical calculations; the second section covers 
mainly experimentation and discussion of measured data versus the cor- 
responding theory. It is recognized that in some cases this division 
of a topic may be a nuisance. However, since part of the material in 
Sec. I has not yet been verified experimentally, and since part of the 
material in Sec. II stands alone and part is related to a theory devel- 
oped in Report 1, this classification was chosen. 

All equations are written in the MKS system of units. However, the 
numerical values of some quantities are given in other metric units for 
convenience as to their size. In most of these cases, only a decimal 
change is required to convert to MKS units. For example, flux values 
are given in CGS maxwells, which can be converted to MKS webers by mul- 
tiplying by 10"" 8 wb/maxwell. 

Section I is divided into four parts. Parts A, B, and C deal with 
leg models and switching parameters and are intended to improve and 
extend our analytical tools in dealing with flux switching in multipath 
cores. In Part A, realistic functions for static 0(F) are calculated. 

In Part B, the concept of a leg model is extended to a va r i ab 1 e - wi dt h 
leg and the ine 1 as t i c - swi t ch i ng parameter p of a tapered leg is calcu- 
lated. In Part C, calculation of average p, which was developed in 


xix 


Report 1 only for A0 * 2 <P r , is extended to the case where A^may vary 
continuously from zero, to 2<p r , Part D contains analytical calculation 
of the f lux- di v ision phenomenon that was described in Report 1 with no 
explanation , Due to the lack of switching models of a leg in a partially- 
switched state, a saturable core (c.g., a shaped MAD), rather than a 
nonsaturable core (e.g., the Laddie) is assumed. Part D is by far the 
most significant part of Sec, I. 

Section II is divided into five parts. In Part A, material prop- 
erties are of concern, and since these properties can be best measured 
by testing a toroidal core, all the experiments in Part A are performed 
on a toroid. In Part B, interest is focused on the geometrical properties 
of a leg and, in particular, on how to divide a multipath core into legs 
and how to determine the leg dimensions in doubtful cases, such as when 
the division into legs depends on where NI is applied, or when the 
switching path does not follow the leg contour. Part C contains test 
procedure and laboratory techniques for measurement of leg parameters, 
which are affected by both the material properties and the geometrical 
properties. In Part D, flux division in a nonsaturable core is studied 
experimentally; experimentation and machine computation of flux division 
in a saturable core are described and compared. Finally, in Part E, the 
unsetting effect, which was described in Report 1, is verified experi- 
mentally. 
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I CALCULATIONS 


The objective of this section is twofold. In Parts A, B, and C, we 
shall extend our analytical tools by calculating the following; a real- 
istic static 0(F) curve, inel ast ic - swi tching parameter p of a tapered 
leg , and average p for partial or full flux switching of a c ons t ant - wi dth 
leg and a tapered leg. In Part D, we shall employ the parabolic model 
for p to analyze flux division in a saturable core. 

A. MODELS FOR STATIC 0(F) CURVE 

In a previous report 1 (hereinafter referred to as Report 1), vari- 
ous models for the static relation between B and H were described. One 
of these models [Report 1, Fig. 4(b)iii, p, 4] approximates this relation 
by straight lines, which are horizontal in saturation and sloping between 
saturation levels. Based on this model, the static 0(F)^ curve of a 
toroidal core was calculated by integrating B from r * r j to r * r o (c/. 
Report 1, pp, 36-39). We would first like to extend this calculation 
to a cons t ant -width leg (c/. Report 1, p. 18). The drawback of this 
model stems from the assumption that B(H) is piecewise linear, when in 
fact B(H) is nonlinear. Our main objective is to calculate 0(F) from a 
more realistic B{H) function, 

1. 0(F) Based on Piecewise Linearized B(H) 

The piecewise linearized B(H) curve between saturation levels 
[Fig. 4 (b ) i i i , Report 1, p. 4] for // M1X > H > H min is described by 

B - ~B r + — ( H - H mln ) . (1) 

max _ i n 

♦ 

References are listed on p. 153. 

^ The term "<f>~F ” in Report 1 is replaced by the term **<t>(F) n in order to emphasixe that 0 is a function 
of F. Similarly, **fl (/f) ” will be used instead of "B-H.” The difference between the terms stat ic " curve ” 
and atatic " loop" is trivial. A stat ic " curve ” is a portion of a stat ic " loop, " and is obtained by 
varying ff (or F) in one direction, so that no hysteresis is recorded. The term curve” in Report 1 

( cf . p. 45) stands for residual (or peak, if so specified) flux. This will be preceded by ** T- microsecond ,** 
e , g. , ”20-M*ec c/>(F) curve.” 
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Consider a leg such as shown in Fig. 19 of Report 1 (p. 18) of constant 
width w, height h f angle Ot and edge lengths and l o , and assume that 
an MMF F is applied along the leg. At a radial distance x from the 
inner edge, 

F 

H - T (2) 

where, following Eq. (26) of Report 1 (p. 18), 


l 


+ ax 



w 


(3) 


Substituting Eqs. (2) and (3) into Eq. (1), we obtain the relation 


B(F,x) 



where 


f m i o 


b r 


'i . 


(4) 


(5) 


and 


pm 




(6) 


The <$>(F) curve is obtained by integrating B{F,x) over the leg cross 
section , i . e . , 


- h j" B(F ,x)dx 
0 


(7) 


The flux state in each of the five segments (Fig. 39, Report 1, 
p. 38) for the case where l o H min < is shown in Fig. 1. The bound- 

aries corresponding to each segment are indicated by dashed lines. In- 
tegrating Eq. (7) over these boundaries, we get the same relations between 
4> and F as expressed by Eqs. (57) through (62) of Report 1 (p. 39), except 
that 2vt . and 27rr are replaced by l. and l , so that f® in anc [ are 

as given in Eqs. (5) and (6), 
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FIG. 1 STATIC FLUX STATES FOR LINEARIZED B(H) MODEL 


2. 0(F) Based on Hyperbolic Models for B(H) 

A hyperbolic model for B(H) in negative saturation was proposed in 
Report 1 (p. 19). According to this model, 



where H a and B ss are empirical material parameters [replacing Parameters a 
and B a in Report 1 (p. 19)]. This model describes the static B(H) for 

H < II where H , is the threshold value of H above which inelastic 

— t h ’ t h 

switching takes place. 
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For H ff th , the following hyperbolic model is a good approximation. 


B = B 

U 


-B, + (B. . + B ) 


H ~ H 

9 

II - H 


f J -o H 


(ID 


where //^ and // r are material parameters. 

The plots of these hyperbolic models for the material B(H) are shown 
in Fig. 2. The intersection point, Point P t represents the transition 
from the region of elastic switching to the region of inelastic switching; 
hence, H at Point P is the threshold field // th . By equating Eqs, (10) 
and (11), we find that 



( 12 ) 



FIG. 2 HYPERBOLIC MODELS FOR STATIC B(H) 
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where 


II 


= "„ + "„ + "n + 7T 

a q n g 


+ ", - "J 


(13) 


Note that the plots in Fig. 2 describe only the material contribution to 
B. By equating B u ~ /J, 0 H in Eq. (11) to zero and B r , we obtain expressions 
for the coercive field // and for ff : 

c r 


and 


// = " + — (". - " ) 
c <7 g 7 n 


H 


(B„ + B )H ~ 28 // 


B ~ B 

s s r 


(14) 


(15) 


The hyperbolic models for B(H) [Eqs . (10) and (11)] are used in 
Appendix A to calculate a static curve of a const ant -width leg. 

The resulting expressions, Eqs. (A-3), (A-6), and (A-7), describe static 
<p(F) in three regions of F : F < H t ^l H t h Z i < F < H t h Z o , and H t h l o < F . 
For low- threshold ferrite material, if | F | is lower than, say 1% of 
B 1 .//x 0 , some terms in these expressions are negligible. These negli- 
gible terms include the so-called air flux and, in the region H th l i £ 

F < H th l o , the elastic flux change in the leg section that has not 
switched ine 1 ast ical ly . By ignoring these terms, Eqs. (A-3), (A-6), and 

(A-7) are reduced to the following expressions: 


for F < H .1 . , 

— t n i 


<P 


hw(B s) ~ B r ) 

(T - 1 , )// 

o la 



- B r hw 


(16) 


! 
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for H t .1 . < F < H. . I 

t h i — — t h o 


4> 


+B r )H q \ 

(I - L)H h 


F /ll 

l.+Fl 

* \H H 


'th 


« i 


In 


H 


F 1 


H 


th y 


'“Vi 


B r hw 


(17) 


and for H r . I < F , 

t n o — ’ 


<P 


M8„ +B r )H q 

(I -LK 




i n 

H H ) ln \F ~ H I., 

> n qf * n i / 


B.hw 


(18) 


Equations (16), (17), and (18) hold for a 1 arge practical range of 
F. We have assumed this range to be ±0.01 so that B is 

less than 1% of B, , . For cores to be described in Sec. II, this amounts 
to about 50 times Beyond this range, the air flux becomes in- 

creasingly significant. 

For |F| >> H th l Q , Eqs. (A-3) and (A-7) may be approximated by 


0 = —<b + cb 

^ 3 % r nr 


(19) 


and 


0 = 0 + 0 . 

33 r a ir 


( 20 ) 


where 


4>. . - B.hw 


( 21 ) 


and 


hw F- 0 (l, 

' F — — u r 1 

o i * 


( 22 ) 


The term </ air stems from the flux density fj. Q H in the cross-sectional 
area hw of the leg. 
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Experimental verification for Eqs . (16), (17), and (18) will be 

given in Sec. II-C-3. 

3. Summary 

Two empirical hyperbolic functions have been proposed for describing 
the static B(H) curve, requiring five material parameters: B r , , 

H , and H . The first function, Eq. (10), covers the region of negative 

n q 

saturation, and includes the parameters B sf , B r , and H a . The second 
function, Eq. (11), covers the nonsaturation region as well as the re- 
gion of positive saturation, and it includes the parameters B,,, B r , H n , and 

H . Based on these functions, the static <p(F) curve was calculated by 
9 

integrating B over the leg cross section. This curve is composed of 
three portions that are described by Eqs. (16), (17), and (18). Calcu- 
lated and measured <fi(F) curves will be compared in Sec. II-C-3. 

B. TAPERED LEG 

Switching models for a cons t ant - width leg were developed in Ref. 1. 

In some commercially available multipath cores, leg widths are not con- 
stant. For example, the widths of the inner legs (Leg 2 and Leg 3) in 
a three-hole MAD [Report 1, Fig. 1(b), (p. 2)] are variable. Another 

example is the Transfluxor, which is an unshaped three-leg core, similar 
to the core in Fig. 123 of Report 1 (p. 126). In the Transfluxor, all 
three legs vary in width. The question is, what are the switching param- 
eters of a variable- width leg, and in particular, how do these parameters 
differ from those of a c ons t an t - wid t h leg? The answer to this quest ion is 
not an easy one. The objective in this section is to discuss the factors 
that are involved in this problem and to propose a partial solution for 
a simple case, namely, to calculate the ine last ic- switching coefficient 
p of a tapered leg. Before we examine a var i able -width leg, it is ap- 
propriate to discuss further the limitations of a leg model in general, 
i.e. , limitations that apply to both const ant -width and var i ab 1 e- wi d th legs. 

1. Leg-Model Limitations 

The difficulty in defining the boundaries of relatively short legs 
was regarded in Report 1 (p. 3) as a limitation on the leg model. The 

leg model has additional limitations that will be discussed as follows: 
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(1) A given F is applied across a leg of edge lengths l ^ and 

l 0 . Let us denote the variables associated with the short 
and long edge by subscripts i and o, respectively. Since 
l ^ ^ H o and, hence, A B i >A B q , This means that 

the flux switching is not uniformly distributed. A leg 
model should account for such a nonuniformity by assuming 
a model for B versus H t and integrating the variation in B 
across the leg width. In Report 1, this was accomplished 
for the static characteristics (pp. 36-39) and for the 
e 1 ast ic - switching model (Report 1, Appendix B, pp. 187- 188). 
This was also done in Part A of this report in calculating 
a static 0(F) curve, based on hyperbolic models for B(H). 

On the other hand, variation in B across the leg is not 
accounted for by the inelast ic - switching parabolic model 


- P(F - F 0 ) 



(23) 


In order to account for such a variation, the parabolic 
model for inelastic switching should be first written in 
terms of B and H t rather than 0 and F, i.e, , 


B * £(H ~ H 0 ) 




(H - H Q ) 


(24) 


where, following Eq. (49) in Report 1 (p, 30), 


4 


p 



(25) 


In Report 1, the parabolic model for p was written as p * p^[l - (0#> r J 2 ] . As pointed out in Report 1 
(p. 28), this model implies the absurdity of infinite switching time for <f> to reach a peak value, unless 
the initial flux level is assumed to be -0.g. slightly above , where <f> rQ Experimental obser- 

vation also indicates that the initial value of p is slightly above zero. In the present report, we 
modify the parabolic model intop = P p [l - (4>/4> s J 2 ] , where 4> $ > 4> r (note that <f> f <f> ). The value of 
4> s is so determined as to yield good fit within a range of excess MMF of application, say up to SFg. In 
this way, the initial flux level, denoted by could be assumed to be -0 r . This assumption is sup- 
ported by the appreciable size of domains of reversed magnetization^ (positive domains) at the static 
state of 4> * -0 r . Additional refinement is based on our assumption that inelastic switching is preceded 
by elastic switching (e/. Report 1, p. 14), and therefore the initial flux level <f> ^ might be above -0 , 
especially when s high F with a short rise time is applied. 
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Equation (24) should then be integrated across the leg 
width. Calculation of this sort for a toroidal core can be 
found in Ref. 2, but the results are too complicated to be 
practical. With the use of a modern computer-programming 
language, such as ALGOL or FORTRAN, this complicated model 
may become useful. 

(2) Due to the nonlinearity associated with the parabolic 
model, Eq. (24), the nonuniformity in B distribution across 
the leg cross section becomes worse prior to the time when 
£.(t) passes through zero. To understand this, suppose 
that the leg is initially at negative saturation, and a 
constant MMF F is applied across the leg. As explained in 
(1) above, A 6. >A B o , and, hence, B i > B q , Substitution 

of the condition £. > B q into Eq. (24) results in B i > J5 o , 
which, in turn, causes Afi. to exceedAB Q even further. 

This nonlinear process continues until reaches peak 

(when B t - 0). Beyond this point, £L decreases; the peak 
of the factor [l “ (B/B # ) 2 ], i.e., unity, moves radially 
outwards as time progresses. This peak reaches the outer 
edge when [l - (B^/B^) 2 ] is less than unity. By this time, 
the parabolic model has tended to equalize the distribution 
of B across the leg cross section, as compared with the 
time when B i * 0. Thus, if a leg is initially at negative 
saturation and a positive MMF switches the leg to a par- 
tially switched state where B i is around zero, the distribu- 
tion of B across the leg cross section is nonuniform (B^ > B 0 ) 
for two reasons that are embedded in Eq. (24). The first 
reason is the factor ( H “ H 0 ) ; the second reason is the 
factor [1 - (B/B s ) 2 ]. The second reason causes B to become 
more uniform at a later time during switching. 

(3) In describing a leg model, it is assumed that the distri- 
bution of B appearing at a leg junction is uniform. The 
validity of this assumption depends on the path of flux 
closure in the rest of the core. If, for example, the leg 
under discussion is shunted by a short and wide leg, this 
assumption may be misleading. In this case, a higher change 
in B would occur in the inner edge of the shunting leg, so 
that the distribution of B at the leg junctions would be 
highly nonuniform. 

(4) The leg model ignores flux leakage due to pole distribution 
on the surface of the leg. We know, however, that an asym- 
metric position of MMF drive will cause a pole distribution 
on the leg surface (see Report 1, pp . 128-138), and hence 
flux leakage will result. As concluded in Report 1, cal- 
culation of the leakage flux is too complex a problem to be 
practically worthwhile. If a leg is very long and thin, the 
flux leakage may make the leg model invalid. 
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are 


Additional limitations on the model for a v ar i ab 1 e - wi d t h leg 
discussed nex t . 


2, Additional Model Limitations of a Var iable -Width Leg 

Consider a variable-width leg of constant thickness h. The maximum 
residual flux level of the whole leg is equal to B f A min , where A min is 
the constricting cross-sectional area. Throughout the leg, except at 
the constricting region, the leg is partially switched, i.e. t 1 0 1 < 0^ . 
Depending on its history, this state may be soft (see Report 1, 
pp. 25-26 and pp. 48-50). The higher is the previously applied MMF, the 
larger the number of domains of reversed magnetization is, and the softer 
the state is. Unless previously driven slowly, a vari ab 1 e - wi dt h leg is 
most likely to be in a soft state; here is where difficulties lie. These 
difficulties stem from our lack of sufficient information about this 
complex state, and are listed as follows: 

(1) We have discussed previously the general limitations of a 
leg model from the point of view of nonuniform distribution 
of B , assuming that the leg is initially in a hard (satu- 
rated) state. Some portions of a variable-width leg, on 
the other hand, are left in a soft state even after a CLEAR 
pulse. These portions are likely to be found along the 
longer switching paths of the leg. However, the soft state 
is characterized by a lower threshold; hence a lower // is 
needed for a given change in B. This factor, therefore, 
opposes the original effect of the path-length ratio in 
making A >&B o . The extent to which this factor makes B 
more uniform is still to be found. 

(2) Static characteristics of a var i ab le- width leg depend di- 
rectly on soft threshold. We need to know more about soft 
threshold before we can calculate the static 0(F) loop of a 
variable-width leg. 

(3) Elastic-switching models have been proposed in Report ] 

(pp. 18-23) for a leg in saturation (i.e., \<p\ > <p ), but 
not for a pa r t i a 1 1 y - s w i t c bed leg. Here, again, the fact 
that a variable-width leg is unsaturated (unless biased by 
very high MMF, in which case surface poles should be ac- 
counted for) prevents us from calculating elastic switching. 

(4) Consider the parabolic model for inelastic switching in a 
constant-width leg, Eq. (23). In looking for a similar 
model for a var i ab le- wid t h leg, we encounter the difficulty 
of the dependence of F Q on the softness of the state. 



Thus, the only switching parameter of a variab le-width leg left to 
be calculated at this time is the inelastic switching coefficient p . In 
order to calculate p of a va ri ab 1 e - wi d t h leg, we must first determine the 
equi-F lines and flow lines. The general solution to this problem is very 
complex ; however, solutions in some special cases are feasible. As an 
introduction to the problem, let us first consider a var i ab 1 e - wi d t h leg 
made of linear material. 

3. Var i ab le - W i dth Leg of Linear Material 

For linear material, the relation 

B - pH (26) 

holds, where the permeability p is constant. The problem of determining 
the distribution of B and H is a classical two-dimensional field problem 
that may be solved by graphical, numerical, and (sometimes) analytical 
methods. An excellent coverage of this topic can be found in several 
textbooks . 4,5 

In linear material, reluctance R is defined as 

R - — . (27) 

0 

It should be emphasized that the term reluctance is entirely different 
from p or p“ 1 . The units of reluctance are (turn 2 /ohm second), whereas 
the units of p are (ohm/turn 2 ). 

Since our prime interest in this work is focused on nonlinear square- 
loop material, rather than on linear material, it is felt that a detailed 
calculation of R of a vari able- width leg should be deferred to an appendix. 
A general expression for R of a leg, whose boundaries can be described by 
orthogonal curvilinear coordinates, is derived in Appendix B. Examples 
of typical problems are also included in this appendix. 

4. p of a Tapered Leg 

We now come back to the original problem, and assume a square-loop 
ferromagnetic material. Calculating p for a variab le-width leg of gen- 
eral shape is much more complex than calculating R, because of the 
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dependence of p on B [Eq. (24)]. Suppose that the effect of the flux- 
closure path on the distribution of B at the leg junctions is negligible. 
If H were invariant along each equi-F line, then, following Eq . (24), B 

would be uniformly distributed across the leg width, and the calculation 
of p would become much simpler. This condition is satisfied in the case 
of a tapered leg. Calculation of p for a leg of this shape is our 
next goal. 





*A-3696-t«4 


FIG. 3 A TAPERED LEG 

The top view of a tapered leg of constant height h is shown in 
Fig. 3. The leg edges, A-B and C-D , are orthogonal to the leg ends, 

4-C and B-D. The geometry of the leg may be expressed in terms of three 
parameters: angle /3, small radius r flf and large radius (the leg 

length is equal to r ^ ~ r a ). shall assume that the band near the 

constricting end 4-C is not driven far into saturation. 

We now divide the leg into curvilinear squares. The equi-F lines 
fall on concentric circles and the flow lines are radial. The radial 
lines meet at the center of the circles, Point 0. Each curvilinear 
square is the top view of a leg element of thickness h, length dr, width 
rd6 , and cross-sectional area hrdd . Multiplying Eq. (24) by hrd6, and 
substituting // = dF/dr , we get 
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(28) 


Bhrdd 


hrdd 
ip dr 



(dF ~ H 0 dr) 


The term ( dF ~ H^dr) is the excess MMF applied to the leg element, and 
the term Bhrdd is the resulting 0 in the element. Following the rela- 
tion 0 = PF ex > P °f the leg element, denoted byA/O, is 


A P 


hrdd / 




(29) 


The inel as t i c - swi t ch ing coefficient of a larger arc-shaped leg element, 
obtained by moving the previous element at a constant value of r from 
8 ■ 0 to 6 * /3, is 


Ip * ftp 

0 



(30) 


Note that B and B s of Eq. (29) have been multiplied by the cross-sectional 
area, h/3r . Since the whole leg is composed of a series of such elements, 
the over-all p of the tapered leg is 



The value of </> # of the leg is determined by the constricting area, i.e 

<*>, - B,hr ai 8 . (32) 

Substituting Eq. (32) into Eq. (31), we obtain 
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Pt 



(33) 


We conclude that for a given <p $t p % is proportional to the angle j3 . We 
also conclude that p t is symmetrical with respect to 0, reaching a peak, 
P tp * at 0 s 0. The expression for p t of Eq. (33) can be written in a 
more general form as a product of p (p and some function 7) t that describes 
the variation of p with (0/0^), Thus, 

p t m Ptp • p, < 34 > 

where 

Ptp 

and 

2 In 


i - (±v 






Let us now compare these results with p of a cons t an t - wi d th leg that 
has the same value of and comparable dimensions: the same height h , 

width w » f3r a (the constricting width of the tapered leg), and length 
l = r b r a , For this cons t ant -wi dth leg, following the relations 
P p * and P = “ (0/0 5 ) 2 ], we obtain 

P m P p ' V (37) 
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where 


and 



( 38 ) 



(39) 


Equations (37), (38), and (39) could also have been derived directly 

from Eq. (33). Referring to Fig. 3, let us imagine that both r a and r fc 
increase indefinitely, while at the same time r fc _ r o = l and r o /3 = w are 
kept constant. Substitution of the resulting relation ( r b /r a ) » 1 + ( l/w )/3 
into Eq. (33), and letting - 0, we obtain (using ^’Hospital’s rule) 

Eqs . (37) through (39). Note that the tapered leg has become rectangular. 

In comparing Eq. (35) to Eq. (38), and Eq. (36) to Eq. (39), since 
> r o , we conclude that (/5 (p /p p ) > 1 and that l 7 ^/ 7 ?) thus, ( P t /p ) > ) 

for any value of 4>/4> s • A more detailed comparison of these parameters is 
shown in Figs. 4, 5, and 6. In Fig. 4, calculated p tp and P p (both in units 



FIG. 4 EFFECT OF r b /r Q ON AND Pp 
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FIG. 5 EFFECT OF r b A a ON r, f AND tjvs. <fi/<p s 



-1.0 -0.5 0 0.5 1.0 



FIG. 6 77/7 VS. 4>/<t> s WITH r b /r a AS A PARAMETER 


16 





of i hP) and the ratio P tp /P p are plotted versus r 6 /r fl . In Fig. 5, cal- 
culated T) (for various values of r fc /r fl ) and p are plotted versus ( P/ ( P S - 
In Fig. 6, the ratio r)Ji) is plotted versus <p/4> s , with r 6 /r Q as a parameter. 

5. Summary 

The inelastic coefficient p of a leg, which is a proportionality 
factor in the expression for 0(F), is a product of two factors. The first 
factor, p t is a function of the material properties and the geometry of 
the leg. The second factor, 7} f is a function of the normalized flux level 
and, if the cross-sectional area of the leg is not constant, also a func- 
tion of the leg geometry. If the constricting end of a tapered leg re- 
mains unchanged as the other end is enlarged, then flux switching becomes 
easier and p % increases. This increase in p t is manifested by the in- 
crease in both p tp and rj t . 

C. AVERAGE VALUES OF p 
1. Introduction 

Average values of p were calculated in Report 1 (pp. 28-31). The 
only reason for using an average, constant value of p is to simplify cal- 
culations where p is involved. As indicated in Report 1 (p. 54), the use 
of average p instead of p(<p) may be a good approximation only whenA0 is 
calculated, not 0(t) or 0(t), because calculation of A0 involves time 
integration of <p{t). 

Two kinds of average p were distinguished in Report 1: t ime - aver aged 

p t p [Eq. (C-7), p. 192], and flux-averaged p, p^ ay [Eq. (C-ll), p. 192]. 
From here on we shall retain the symbol p , but replace P^ By hy P av for 
the sake of brevity. Our first goal will be to examine the meanings of p 
and P , and consider their usefulness. 

In Report 1, we calculated p and p ayf assuming a complete switching. 

In many applications, flux switching is only partial, i.e., the flux varies 
from a CLEAR state to some final value <fc ^ between. ”0 r and +<fi s . In these 
cases, an average value of p is not constant, but rather a function of 0y* 
Our second objective is to calculate the average p for a leg whose flux 
level varies from ~<p r to (p^ , assuming a parabolic switching model for the 
material [Eq. (24)]. We shall perform this calculation for a constant- 
width leg and a tapered leg. 


17 



2. p AND ^ ay 

Suppose that an MMF pulse F of an arbitrary shape is applied along 
a leg in the positive direction. The resulting inelastic flux change 
isAc£ * <t> f - <p Qt where <p Q is the initial flux and <j> is the final flux. 
The inelastic-switching coefficient p varies during the switching of <f>. 
In an attempt to average p over we consider two kinds of averaging 

process: (1) averaging p itself, and (2) averaging the inverse of p. 

The flux- averaged p is defined as 


P 


a v 



pd<t> 


The flux- averaged (1/p) is defined as 


(40) 



(41) 


We now want to determine which of these average quantities may be applied 
in calculation of flux switching. Following the basic definition of p 
[ cf . Eq. (23)], <p ■ pF ex , where f ex (* F - F 0 ) is the excess MMF over the 
dynamic threshold. Based on this relation, 


(■“)# a f F ' X dt (42) 

J ' p ' o e x 

0 

where r is the switching time. Note that if the switching is time-limited 

(as opposed to flux- limited) , then r is the time during which F > F n . 

Because of its dimensions, the term J F dt is called "excess charge- turns 

0 e * 

or in short, “excess charge.” Combining Eqs. (41) and (42), 



( 43 ) 


18 



Thus, in order to compute A0 that results from a drive of a given amount 
of excess charge, we have to use the flux- averaged value of (1 / p) , 

Eq. (41). On the other hand, the flux- averaged p [Eq. (40)] does not 
relate F fx to the resulting change in 0. We, therefore, conclude that 
(1 /p) (and not P ay ) is the useful average value for calculation of p. 

In Report 1 (Appendix C) , p was calculated for switching from 
0 o « -<p r to = +0 r due to a constant F. Examination of Eq. (C-8) 
[Report 1, p. 192] leads to a more general definition of p, i.e.. 


P 



(44) 


A comparison between Eq. (43) and Eq. (44) reveals that 1 Ip * (1/P>„. 
Following Eq. (41), we conclude that 


P = 



(45) 


In Report 1, p was called 11 time- averaged p, ” yet Eq. (45) indicates that 
p is the inverse of the flux- aver aged (1/^). We observe, however, that 
Eq. (44) may be written as 


A0 « PF ex r 


(46) 


where F is the time-averaged F . For this reason we shall continue to 

ex D e x 

call p " time-averaged p r ” although its evaluation will be based on Eq. (45). 

We have concluded that p, not p ay , is the useful parameter. However, 
if p does not vary too much during switching, then to a rough approxima- 
tion, ] /p ** (l/p) #y . Under this condition, p may be approximated by 
p #y (c/. Report 1, p. 29). 


3. p of Constant -W idth and Tapered Legs 

We shall now use Eq. (45) in order to calculate p of a cons t an t - w id th 
leg and p of a tapered leg (to be denoted by p f ) . first, let us 
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rewrite Eq. (45) in terms of normalized flux u, where u = < p/<p g . In ad- 
dition, initial value of u, (= <p^/<p s ) , and final value of u, Uf 

(a <j>f/4> s ), will be spelled out. Thus, Eq. (45) is rewritten as 


P 



(47) 


a. CALCULATION OF p 

The parabolic model for p of a cons t ant - width leg is [ cf . 

Eq. (23)] 


p = p p (1 - u 2 ) . (48) 

Substitution of Eq. (48) into Eq. (47) gives 

u f ’ u 0 

p 9 p p — — : 3 — . ( 49 ) 

tanh "* 1 \if - tanh "" 1 u Q 

Let us examine Eq. (49). For 0 = u f = u r (= <fi r /<p s ), Eq. (49) is re- 
duced to p * P p u r / tanh 1 u r , which, as expected, is identical with p for 

similar switching conditions [Report 1, Eq. (43), p, 28]. For u 0 =* 

Uf = "tt r (no inelastic switching), Eq. (49) is reduced to p * ( 1 ~ u 2 ) 

by using L’Hospital's rule; this result agrees with the parabolic model 
for p [Eq. (48)], in which u = ™u r . We note that p -• 0 as Uf 1, This 

can be explained by referring to Fig. 30, Report 1, and to Eq. (44). As 
Uf — 1 , T and sinceA0 is finite, p -* 0. 

In most cases, a leg is initially in a CLEAR state, so that 
u Q = “u f . Under this condition, Eq . (49) is modified to the following. 

“/ + U r 

P = P p 3 . (49a) 

tanh 1 u, + tanh -1 u 
f r 
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Note that p(u^) reaches a maximum. Assuming that u f (= < P r / ( P s ) = 0.9, 

p/p vs, u, (= <p,/<p ) is calculated from Eq. (49a) and plotted inFig, 7. 
P 7 7 * ^ 

A maximum value of ( p/p ) = 0.695 is reached at *= 0.55. Since 

u Q = _u r » there is no physical meaning to the plot where u ^ 
hence, the region “1 £ <_ “0.9 is marked by dotted line. 



BA-3696-IS5 

FIG. 7 p/p p vs. 0 f /0 s FOR A CONSTANT-WIDTH LEG 

The value of u r in Eq. (49a) must be known a priori in order 
to calculate p . Unfortunately, u r may vary (around ~0.9) among dif- 
ferent materials. The effect of u r on p is shown in Fig. 8, where p/p p 
is plotted vs, u ^ with u r as a parameter. As u r increases from 0.80 
to 0.98, the maximum value of {p/p p ) decreases from 0.79 to 0.53 in the 
range 0.48 < u ^ <_ 0.7. Where dotted lines are drawn, the calculated 
curves have no physical meaning. 
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A 


RB-5696-16* 


FIG. 8 PLOTS OF p/p p vs. WITH <p/<p s 

AS A PARAMETER (Constant-Width Leg) 


6. CALCULATION OF p 

Substitution of Eq . (33) into Eq. (47) (c/. Appendix C) results 
in the following expression, assuming that = ~u r . 


p, - 


u, + u 

7 r 


C+Ur 1“U, C+Ii I“U 

(c + U/ ) '(1 - U/ ) f (c + u p ) r (l - U r ) ' 


(50) 


in 


C“U, 1 +U , C - U 1+U 

-(c - u f ) Hl + U f ) f (c-u r ) r (l + U p ) ' 


where 


c 



As Uf /O f -* 0 (due to infinite switching time). For the case in which 

~u 0 = Uf = u r , Eq. (50) is reduced to 
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P t - 2C//3 


(50a) 


In 


. C ^ u 1 “ U — 1 

(c + u.) r (l - u ) r 


c - u 1+u 

(c - o r ) r (l + u r ) r 


and for u r =* 0,9 it is further reduced to 


Pt 


1.8 i p h/3 


In 


0.2344 


(c + 0 . 9 ) e + 0 ' 9 "1 
(c - 0. 9) c_0 ' 9 J 


(50b) 


Let us check Eq. (50) by letting r a “* 00 and r fc — 00 in such a way that 
l (= r. - r ) and w (= r 0) are kept constant. Substituting c 3 1 + ( l/w)@ 
into Eq. (50) and letting /3 -* 0 , the use of L'Hospital's rule results in 
p being equal to p of a constant-width leg, Eq. (49), as expected. 

Assuming that u r ’ 0.9, P t /P tp versus <P f /4> s is calculated from 
Eq. (50) and plotted in Fig. 9 for different values of c (= r fc /r a ). As 
c increases, Pt^Ptp i ncreas es due to the increase of with c (c/. in 

Fig. 5). 



FIG. 9 p t /p vs. cj&j/c p s WITH r b /r Q AS A PARAMETER (Tapered Leg) 
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4. Summary 


Compared with f 1 ux- averaged p (denoted by P ) t t i me - a ve r aged p 
(denoted by p) is essentially the inverse of fl ux- averaged (1 Ip). In 

T 

calculation of A 0 due to an excess charge J F dt , p (and not p ) is 

0 c x a v 

the useful parameter. If p does not change by a relatively large amount 
during switching, p may be approximated by , 

Calculated plots of p/p p vs. 0^/0^ of a constant- width leg are shown 
in Fig. 8 for various values of 4 > r /4 > 9 * Calculated plots of P t /P tp vs. 
< Pf/ < P t of a tapered leg are shown in Fig. 9 for various values of r^/r fl , 

D. FLUX DIVISION IN A SATURABLE CORE 
1. Introduction 

The problem of flux division 6 (see Report 1, pp. 75-85) is briefly 
as follows: Magnetomotive force NI is applied to one leg in a three-leg 

core, and the resultingA0 in the driven leg is divided into the other 
two legs. Experimental data show that the ratio between the net amounts 
of A0 in the two legs is far from being equal to the inverse of their 
path-length ratio. The question of why this is so and how to calculate 
.such a division of flux in various types of multipath cores presents 
itself. 

Multipath cores may be classified by whether or not all legs in the 
core can be brought to a saturation residual state, i . e . , to 0 = ±0 . 

We shall thereby distinguish between saturable and nonsaturable cores. 

Two extreme cases of such a classification are the Laddie and the shaped 
MAD (see Report ], Fig. 1, p. 2). 

In the Laddie, all legs have equal cross-sectional area; thus, it 
is impossible to have residual states of \cp\ = 0 r in all legs simultane- 
ously. For example, if in two legs of a two-hole Laddie (Report 1, Fig. 1, 
p. 2)0= 0 r , then in the third leg, 0=0. To illustrate another ex- 
ample refer to Fig. 94 of Report 1 (p. 99), in which a three-hole Laddie 
is shown. Although \<p\ = 0 r in Legs 1, 2, 3, and 4, Leg m (comprised of 
the two horizontal middle legs) is in a demagnetized state of 0=0. 

Note that the flux level in Leg m (or the third leg in the previous ex- 
ample) may certainly be other than zero, but then | 0 | / 0^ in the other 

legs. We thus conclude that it is impossible to bring a Laddie (of any 
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number of holes) to a residual state in which all legs are in a CLEAR 
state, M = <p r . A Laddie is, therefore, a nonsaturable core. The 
leg(s) in a partially-switched state of | 0 | / <fi r may have a soft thresh- 
old, which is a complex phenomenon (c/. Report 1, pp. 48-50). Since 
flux division is highly dependent on soft threshold, calculation of flux 
division in a Laddie is a rather difficult task. Experimental data on 
flux division in a two-hole Laddie was reported in Refs. 6 and 7, but no 
attempt was made to explain this phenomenon. 

In a MAD (later, Fig, 42), leg cross sections are so shaped that it 
is possible to bring all legs to a CLEAR state, |0| = <p r . This shaped 
MAD is, therefore, a saturable core. With no legs initially in a soft 
state, the analysis of the flux division is easier than in the case of 
a Laddie. Calculation of flux division in a saturable core was first 
reported in Ref. 8, but was incomplete for the following reasons (cf. 
Report 1, pp. 153-154): First, a square-root parabolic model for p was 

used; although leading to simpler mathematics, this model is less ac- 
curate than the parabolic model for p . Second, it was unjustifiably 
assumed that H in each leg is constant during the switching time. Third, 
rather than being solved for, the values of H in the two parallel legs 
were included as parameters in the final expression for flux division; 
furthermore, the applied MMF to the third (main) leg and the parameters 
of the third leg were completely ignored. Fourth, no experimental veri- 
fication was given. 

Our objective here is to calculate flux division in a three-leg core 
as a function of leg parameters (assuming a parabolic model for p) and 
applied MMF, and to compare this calculation with experimental data. 
Lacking f 1 ux- swi t ch i ng models for a leg initially in a soft state, we 
cannot at the present time analyze flux division in a Laddie. We shall, 
therefore, narrow our scope in this report to calculate flux division in 
a saturable core whose legs are initially in a CLEAR state. However, in 
order to be more general, we shall assume that the output leg is loaded 
by resistance R l across N L turns; flux division in an unloaded core would 
later be determined as a special case of (/V^//? L ) * 0. 

2. Assumptions 

The analysis of the flux division is based on the following assump- 
t i ons (see Fig . 10): 
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FIG. 10 FLUX DIVISION IN A SATURABLE THREE-LEG CORE 


(1) The core is saturable. Thickness h is common to all 
legs, and w 3 = w 4 = w m / 2; hence 

*r. * * », • <5D 

(2) Resistance /? t loads Leg 4 across turns. 

(3) The core is initially in a CLEAR state, t.e., 

= ~ 2(p r - ^3 = ' 

(4) The SET current pulse, applied to N turns on Leg m, is 

rectangular and of amplitude/ andduration T so adjusted 
that the resulting flux change in Leg m , , is equal to 

2 0 . [This assumption of time-limited switching could have 
been eliminated by assuming that NI was applied to Leg 3 
of a five-leg core, as in Fig. 102(a) of Report 1 (p. 107). 
Such an alternative assumption would make the calculations 
more complicated without helping to understand the mecha- 
nism of flux division. We prefer, therefore, to maintain 
the assumption that NI is applied to Leg m and is adjusted 
both in magnitude and duration to yield A0^ = 2<£ r . ] 
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(5) Since flux division is defined in terms of residual 

flux, we shall neglect the difference between A <p e dur ing 
the rise of I and during the fall of I, and consider 

only inelastic switching. 

3. Calculation of Flux -Div is ion Ratio 

The method of calculating the flux division (A0 3 /A0 4 ) is outlined 
as follows. An electrical-circuit analogue is drawn, and node and loop 
equations are written. The parabolic model for p, p - p^ [1 — (cp/4 > 9 ) 2 ] , 
is employed in the equations for 0 of each leg. In this way, 0 4 gets 
to be smaller than 0 3 not only because l 3 < i 4 (and hence // 4 < fl 3 ), hut 
also because p 3 > p 4 during most of the switching time as a result of 
0 3 being higher than < p^. The flux division is expressed in terms of the 
final flux level of Leg 3, 0 3 f. The remainder of this calculation is 
devoted to solution for 0 3 ^. and to expression of A0 3 /A0 4 in terms of NI 
and the leg parameters. 

Based on Assumption (5), only the inelastic parameters Fq and p of 
each leg are considered. The general electrical-circuit analogue for the 
loaded core in Fig. 10 (c/. Fig. 37, Report 1, p. 35) then becomes simpler, 
because all the “capacitors” are omitted. Since only the SET time is of 
interest, the e 1 ec t r i c a 1 - c i r cu i t analogue is further simplified by omitting 
redundant components (c/. Report 1, Fig. 37, p. 35) as shown in Fig. 11. 
Since in an e 1 ec t r i c a 1 - c i re u i t analogue MMF is treated as “voltage” and 
0 is treated as “current,” the SET MMF, /VI, appears as a “voltage source” 



FIG. 1 1 ELECTRICAL-CIRCUIT ANALOGUE 
FOR A DRIVEN LOADED CORE 
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in series with Leg m, and the back MMF, appears as a “voltage drop" 

across a 11 resistance n N^/R L in series with Leg 4. A special case is that 
of no load on Leg 4, t.e., N^/R L * 0; this case can be obtained by letting 
either R^ 9 & (open circuit) or B 0 (no coupling). 

One node equation and two loop equations are written by inspection 
of the electrical-circuit analogue, Fig. 11, as follows. 


4 > m • 0 3 + <*> 4 


(52) 


and 


NI = F + F, 

m 3 

"l • 

3 Fa " 7 


(53) 


(54) 


Using a parabolic model for p [cf. Eq. (23)], 


and 



(55) 


(56) 


(57) 


Equations (52) through (57) contain six unknowns: , <fi 3 , <p F # , F 3 , 

and F Although these unknowns can be solved for, there is no need for 
these solutions in calculating the flux division, as will be shown next. 

Let flux-division ratio be defined as 


A03 

A 0 4 


(58) 
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where A0 3 andA0 4 are the inelastic changes in flux in Legs 3 and 4, 
respectively. Since both of these legs are initially at 0 * “0 r , 

A0 3 - 0 3/ + 4> r < 59 > 

and 

A0 4 - <fi tf + 0 r (60) 

where subscript / denotes final flux level. Following Eq. (52) and 
Assumpt ion (4 ) f 

A0 3 +A0 4 “ A0 k = 2 4> r . (61) 


From the substitution of Eqs . (59) and (60) into Eq. (61 ) t we con- 

clude that 


= ~ <p *f 


(62) 


Typical variations of 0 versus F of the three legs during the SET time 
are plotted in Fig. 12, superimposed on the corresponding static 0(F) 
loops. These plots illustrate why 0 3/ » as stated in Eq. (62). 

Substituting Eqs. (59), (60), and (62) into Eq. (58), the flux-division 
ratio is found to be 


D 


0 


1 + 


3 / 


4>, 


0 


1 - 


3 / 


0 r 


(63) 


Our next objective is to calculate 0 3 ^. 

During the switching time T , 0 3 changes from ~ < P T to 4 > ^j ( c / • Eig. 12) 
Following Eq. (56), 


3 / 


d0, 


A \ 2 


1 - 





3 p 


i (E 3 - F 03 )dt 


(64) 
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t*. 



FIG. 12 VARIATION OF q> AND F OF EACH LEG IN 
FLUX-DIVISION EXPERIMENT 
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After integration, we get 


(f> 


- 1 

t y - + 1 


in which 


7 = 2 


(F 3 - F' 0 ,)dt ~ tanh"^ — 


Substitution of Eq. (65) into Eq. (63) gives 


1 - Se r 


in which 


<t> r 
1 + 4>. 


Calculation of y is performed in Appendix D. It is shown there that 


y 3 p ~ q In S 


where 


"i\ *, t 

(NI - F 0m )(p Sp + P 4o ) ” ~ 


04 ^ 4 p 1 0 y 3 p 


w - F 0m )ip 3P - P 4p ) + F 0 ^ 4p ^1 + 2 — y F 03 p 3p ^1 + 2 — J 
a * (71) 

w - + />„> - f 04 p„ - f 0 ^„ 
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Note that the expressions for p and q have the same denominator. The 
value of q is independent of iV^/7?^ an d can be determined directly from 
Eq. (71), but the expression for p [Eq. (70)] includes the st i 1 1 - unknown 
term . By combining Eqs . (65) and (69), we obtain 


til 


t anh 



q In 

2 


(72) 


After q is determined, substituting Eq. (72) into Eq. (70) results in a 
transcendental equation in p . A numerical solution for p can be obtained 
without any difficulty. 

Substitution of Eq. (69) into Eq. (67) results in the final expres- 
sion for the f lux- di vi s i on ratio: 


D 


ef - 5 9 * 1 

5 ^ - 5 e p 


(73) 


in which 5 is expressed in Eq. (68), p is expressed t ranscendent al ly in 
Eqs. (70) and (72), and q is expressed in Eq. (71), 

The larger the load N*/R^ is, the smallerA<£^ is, and hence D be- 
comes higher. Mathematically, for the same values of NI and core param- 
eters, if NI/R l increases, then p [Eq. (70)] increases as well; hence, 
e p increases. An increase in e p causes D [Eq. (73)] to increase because 
of an increase in the numerator and a decrease in the denominator. 


4. Flux -Div is ion Ratio Under Special Cases 

The general expression for D [Eq. (73)] will now be reduced under 
special cases. 


a. UNLOADED CORE 

As explained before, if either N L = 0 (no coupling) or R^ a 00 
(open circuit), Leg 4 is unloaded. Substituting {N^/R L ) - 0 into Eq. (70), 
we find that p = 0. Recall that q is independent of /R L . The ex- 
pression for 0 [Eq. (73)] is then reduced to 
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1 - s* 


( 74 ) 


+ 1 


where 8 and q are expressed in Eqs . (68) and (71), respectively. 

6 . D - 00 

When A 0 4 - 0 andA<£ 3 > 0 , D “V 00 . This condition is physically 
realizable, since if is large enough, NI may be so adjusted as 

to cause inelastic switching to occur in Leg 3 only. SinceA<£ 3 +A0 4 * 2 <p r 
andA<^ 4 * 0, thenA0 3 * 2 <p r \ hence, following Eq. (59), 0 3 f a <£ r » » n d 
therefore, p [Eq. (70)] is zero. Substituting D 00 and p = 0 into 
Eq. (73) results in q * 1. By equating q [Eq. (71)] to unity, the value 
of NI for which D 00 i s found to be 


NI 


/)-<*> 


0 M 



+ F 


04 



(75) 


The expression for NI [Eq. (75)] is obviously not a function 
of the loading factor, . This stems from having p, the only param- 

eter where N*/R t has an effect, be equal to zero. The physical inter- 
pretation is rather simple: - IfA0 4 » 0, the presence of a load on Leg 4 
is irrelevant; hence, the magnitude of has no effect. We also 

notice that NIp^^ is not a function of 8 (since q is not a function of 8). 


c . /VI ^ 00 

It was reported 6 that the flux-division ratio in a Laddie core 
approaches an asymptote as NI “* 00 . Let us now calculate the asymptotic 
value of D [Eq. (73)] for a saturable core. 

As NI ~* °°, p and q [Eqs. (70) and (71)] are reduced to the 

f ol lowing : 


Pn I -a 


2 N _l PlpPjp I’K 

\ Pip + PipX^s 



(76) 
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and 


Qni~ CD 


Pip ~ Pip 
Pip + Pip 


(77) 


After i s determined, Eq. (72) is substituted into Eq. (76) and 

p NI ^ © is solved for t ranscendent al ly . The values of p NI ^ca an d <7* j-*cx> are 
then substituted into Eq. (73) to calculate D^j^. 

The larger the load N*/ H L is, the larger is p N ^co [Eq* (76)]; 
hence, following Eq. (73), the larger is The lowest is 

obtained when the core is not loaded. In the no-load case, p 3 0 and 
<7jy/-*co [Eq. (77)] is substituted into Eq. (74) in order to calculate the 
asymptotic value for D NL , i.e., 



(78) 


d . § - l 

Following Eq. (68), if (0 r /0 4 ) “V 0, then 8 -• 1. Before deter- 
mining how D would be affected by the condition 8 1 , let us examine 

the physical meaning of having (< t > r /4 > s ) 0. 

Our analysis has been based on the assumption that initially 
(fi* s ~2<£ r and </> 3 = 0 4 = “0 r . In calculating A0 of each leg, the para- 
bolic model p = p { 1 ~ (4>/<P 9 ) 2 ] was employed. Suppose that <fi r is replaced 
by 4> TeB i where 0 rea is an arbitrary residual flux level, so that initially 
d> =* -2<h and <£, = <t> . - “0 . If this initial state were obtained by 

slow, gentle switching, calculationoi$ in each leg may still follow Eq. (23). 
Under this condition, our analysis still holds, except that 8 [Eq. (68)] 
is modified by replacing <fi by 0 ; hence, 8 1 when (0 /<£ ) 0. 

Modifying Assumption (4), = 2</> res , The physical interpretation of 

the condition {<fi s /4> s ) 0 is, therefore, as follows: All legs switch 
by a very small amount near the demagnetized state of <p - 0, and hence 
the inelastic- switching coefficient p of each leg remains practically 
constant at its peak value, as illustrated in Fig. 13. Mathematically, 
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FIG. 13 p/p p vs. 0 OF EACH LEG WHEN S - 1 

each of the bracketed terms [1 ” (0 a /2$ f ) 2 ], [1 " (0 3 /0 j ) 2 ], and 

[1 “ (0 4 /0 s ) 2 ] in Eqs. (55), (56), and (57), is replaced by unity. 

We now want to derive an expression for D under the condition 
that S “* 1 . The transcendental solution to Eq. (D-21), as (0 r /0 s ) 0, 

yields y — 0. Derivation of D from Eq. (67) and Eq. (D-21) is then 
based on L'Hospital's rule, which gives 


1 im D * 
( 4 > r /< t> s )**0 


a 3 + a 2 + 2a t 



where a t , a 2 , and a 3 are expressed in Eqs. (D-17), (D-18), and (D-19); 

thus , 


W - r o. -F 0 3 )p mp + ( F ot - F 03 )p./[ 1 + (pMf l )} 


Pj^ + 


m - F r 


F 04*., + ^03 


^04^°3p 


(79) 
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It can be shown by simple calculation that Eq. (79) would re- 
sult if constant p was assumed in the original equations [Eqs . (55) 

through (57)], i.e., if p m = p mp , p 3 = p 3p , and p 4 = p ip . 

By equating the denominator of Eq . (79) to zero, the asymptotic 

expression for NI , for which ^ is determined. This expression is 

identical with the expression in Eq. (75), as expected, since the latter 
was shown to be independent of S. 

In the no-load case, (jV^//?^) » 0, and Eq. (79) may be reduced to 


D 


Si 


, 5-1 


/ ^3 p \ ^0n ~ ^ 03 )P mp + (^Q4 ^03^4 p 

\ P ip J W ~ F 0, - F 0 *>P mp + < f 03 ' ' 04 >P„ 


(80) 


When NI °°, Eq. (79) may be reduced to 


0 8-1 = — Cl + (P 4p /V t 2 /fl t )] (81) 

P \ p 

In the no-load case, Eq. (81) may be further reduced to 

® SL ,5-1 j-o o _ (81a) 

^4 p 

In conclusion, following Eq. (79), vs. NI is a hyperbola 

whose two asymptotes are defined by NI * NI [Eq. (75)] and by 
D = D s „ lpjrjM * [Eq. (81)]. 

Let us summarize what we have found so far. D versus NI has 
been calculated in terms of leg parameters. We have seen that for a 
given value of NI , D decreases if the loading factor N^/R l decreases, 
reaching minimum when the core is unloaded, i.e., when N^/R L = 0. The 
asymptotic value of NI , for which D 00 , was found to be independent of 
«V*l or S. lhe asymptotic value of D as NI 00 decreases when N^/R f de- 
creases, reaching minimum when the core is unloaded. This asymptotic 
value of D for N I — 00 also decreases as & increases, reaching minimum 
when S 1, i.e., when in each leg A</; is very small around <£> = 0, and 
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FIG. 14 D vs. Nl IN THREE CASES 

p *= p = constant. All this is summarized graphically in Fig. 14 (assum 
ing that the leg parameters are known) by sketching D versus NI in three 
cases: > 0; = 0; and 8 -♦ 1 in the case of * 0. 

5. Effects of Leg Dimensions 

Up to the present, calculation of flux division has been based on 
the leg parameters F^ m> F 0 ^> P mp * P3p> an< * P 4 p • Our next problem 

is to determine these parameters. Since the three legs are of the same 
material, the leg parameters differ from each other only because of dif- 
ferences in leg dimensions. Determination of the effect of the leg 
dimensions on D, a problem that was first introduced in connection with 
the Laddie, 6 amounts, therefore, to determination of the effect of the 
leg dimensions on the leg parameters. Before we determine this effect, 
let us improve our inelastic switching models in the region of low F. 

a. LINEAR AND NONLINEAR REGIONS OF B VERSUS H 

P 

Material i ne 1 ast ic - switch ing properties have been discussed in 
Part B. Following Eq. (24), when B * 0 , B reaches a peak value of 

K - C p (H - // 0 ) • (82) 
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( 84 ) 


where is the intersection with the H axis and 4* is the slope of the 
tangent to the nonlinear portion of B p versus//, as shown in Fig. 15. Note that 
H q and 4 p may be looked at as a special case of H ^ and in which B p 
versus H is linear. Equating dB p /dH of Eq. (83) to dB p /dH of Eq. (84) 
gives 

4 ; = kp(H - H u 0 ) v - 1 (85) 


Substituting the expression for B p [Eq. (83)] and 4^ [Eq* (85)] into 
Eq. (84) gives 


H 


0 



( 86 ) 


Inelastic switching may also occur for H £ H"^. However, the 
parabolic model for 4 [Eq. (24)] does not hold for this region of slow 
switching. Further investigation is needed in order to determine models 
for such a slow inelastic switching. 


6. GEOMETRICAL RELATIONS 

We now introduce leg dimensions in order to describe inelastic 
switching in terms of <$> p versus F rather than B p versus H. Integration 
of B across a leg cross section could not be expressed by a closed form. 
Because of this mathematical difficulty, we shall ignore the variations 
of H and B across the leg. Replacing // by F/l (for the sake of brevity, 
the symbol l By is replaced by l ) and B p by 4> p /A t Eqs. (82) through (86) 
are transformed into the following: The linear region of 4> p versus F is 

described by 

<4 = P p ir-F 0 ) (87) 

where 

p = 4 — (88) 

p p l 

and 

F 0 = V ■ (89) 
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The nonlinear region of <p p versus F is described by 


K(F ~ F u 0 ) v 


where 


For a given value of F , <p p (F) may be described alternatively by 


where 


~ r' 0 > 


\v(F ~ F " Q ) V ~ 1 


p u 


The relations between the leg parameters can now be derived. 
These relations may later be used to express D as a function of the leg 
dimensions. Let us refer the leg parameters to one of the legs, e.g. , 
to Leg 3. If the legs to be compared operate in their respective linear 
regions of <p p versus F, then, following Eqs. (88) and (89), 


A l, Z, 

p* = £> 2 — 

3p A. I 3p l 

3 at at 


^4^3 ^ 3 

A l ~ y ° 3p l 

/l 3 4 4 4 4 
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(98) 


= F 


03 


and 


*4 

^04 = ^03 i 

1 3 


(99) 


If the legs to be compared operate in the nonlinear regions of <fi p versus 
F, then, following Eqs. (91) and (92), 


and 



(100) 


(101) 


( 102 ) 


(103) 


Before the relations between the leg parameters are substituted 
into the expressions for D , we need to derive the conditions that deter- 
mine whether each leg operates in the linear or nonlinear region. This 
derivation is performed next. 


C. VARIATIONS OF H IN EACH LEG 

Suppose that the cross-sectional areas of Legs m, 3, and 4, as 
well as the lengths of Legs m and 3, are fixed, but the length of Leg 4 
is varied. We shall first assume that l 4 >_ Z 3 , and examine two extreme 
cases: In Case a, l ^ * i 3 ; in Case 6, 1 4 - 00 . We shall assume that 

Leg 4 is unloaded, since loading Leg 4 has essentially the same effect 
on flux division as increasing l of an unloaded core, 8 and thus the 
case of loaded Leg 4 (of finite length) lies between Case a and Case 6. 
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When Z 4 3 Z 3 (Case a), Legs 3 and 4 in parallel appear as a 
single leg whose cross-sectional area is equal to A m . In this case, the 
applied NI would divide uniformly along all legs, and // 3 =* // 4 * H u . 

When 1 4 00 (Case 6), the core is essentially composed of 

Leg m and Leg 3 in series. Substitution of p p = £, p A/l, F 3 HI, and 
F o * H into Eqs. (55) and (56) yields expressions for <P m and <p^ . 

These expressions are equated, since <P m - 0 3 , yielding the following 
relation: 



(104) 


It is evident from Eq. (104) that the relationship between H m and H^ is 
affected by the cross-sectional areas of Legs m and 3, but not by the 
lengths of these legs. Substituting the relations {A^/A m ) = 1/2, 

4> m * “2<p r and 0 3 D "0 r + A0 3 * “0 r + A into Eq. (104), we get: 


K ~ H _o 

", ~ 



<p \ 2 

A0. _ <M 2 

2<P S 4>J 


(105) 


As switching progresses in time, A 0^ varies from zero to 2<fi r . The curve 
( H m “ H 0 )/(H 3 ~ H Q ) versus {&<p*/2<p s ) [Eq. (105)] is plotted in Fig. 16, 
assuming that (4 > r /4 ) s ) s 0.9 [hence (Ac^/20^) 3 0.9 corresponds to 
A<£ # * 20 r ]. Note that the curve is below unity during the whole switch- 
ing time, and that it reaches a peak of ^0.735 when A 0^ = 0.4 <p $ 3 O.4440 r , 
i.e., when Leg m completes about 22.2% of its total switching. 

It was shown previously that in Case a (where Z 4 * Z 3 ), 

H m - H z , and hence (H m ~ H Q ) / (H 3 - H Q ) = 1. The plots of (H m ~ H Q ) / (H ^ ~ H Q ) 
versus (&<P m /24> s ) in the two extreme cases are compared in Fig. 16. Any 
intermediate case of 1 3 <z 4 < ro 

would result in a plot that lies between 
these two plots. Thus , if [ 4 > Z 3 , then H m < H ^ throughout the switching time. 
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'"-CASE a: U*l3 



FIG. 16 (H m - H 0 )/(H 3 - H 0 ) vs. (A0 m /2 0 # ) 

From the way ( H m ~ H 0 ) / (H ^ ~ H 0 ) versus (A0 (| /20 s ) varies in 
Fig. 16 as Z 4 is decreased from 00 to it may be deduced that if l 4 

is decreased below Z 3 , the ratio {H “ // Q ) / ( // 3 ~ // q ) would rise above 
unity. Suppose now that Leg 3 and Leg 4 interchange functions, and we 
plot (H m ~ H Q ) / (H “ // q ) ve r sus A <p* /2<p $ . Then, the curve above unity 
will correspond to Z 3 < l 4 . Applying this observation to our previous 
case in which < Z 4 , we infer that H m > H 4 . In conclusion, for any 
value of NI large enough to switch flux inel as t ic al ly in a three-leg satu 
rable core (A ^ = A^ = A u /2) , if Z 4 > Z 3 , then > H m > H 4 during the 

whole switching time. Since H 3 > H m > H 4 at any instant of switching 
time, we deduce that // 3 > H m > H^, where H indicates a t i me - aver aged 
value of H over the switching time T. 




These inequalities are of interest, but additional information 
about the magnitude of H in each leg relative to H B (Fig. 15) is needed 
in order to determine the region of operation. Since calculation of in- 
stantaneous H in each leg and its effect on D is quite involved, we shall 
simplify the problem by replacing H by H . Calculation of //-, H , and H A 
(assuming that Leg 4 is unloaded) will be done next. 

With no load on Leg 4, 

, T 

{ F > di 

F, ■ F, - — - — • (106) 

T 

Substitution of the expressions for J F~dt [Eq. (D-12)] and T [Eq. ( D- 14)3' 

0 J 

into Eq. (106) gives 


F 


3 



m - 


P + 

r mp 


* [ojPjp + _ [03P3 p 

Pip + Pip 


(107) 


Since F m = NI ~ F 3 , then F m = NI ~ F 3 ; hence, F a can be obtained by sub- 
tracting F 3 [Eq. (107)] from NI . Following the relation H = F/l, we 
calculate H 3 * F^/l^, H i = F i /l i , and H m = F m /l m . Each of these Hv alues 
may be compared with H B , Fig. 15, in order to determine whether operation 
occurs in the linear portion of B p versus H (or 4> p versus F) . It should 
be emphasized that the conclusion drawn by this step is only an approxi- 
mation, because the instantaneous value of H was replaced by its time- 
averaged value, H, 

It was shown previously that if J 3 < 1 4 , then H 3 > H m > H^. 
Therefore, if , then all three legs operate in their respective 

linear portions. This condition is satisfied when NI > NI ot , where NI n 

— B 4 B 4 

is the value of NI for which (F^/l^) 3 H B . Under this condition, the 
relations between the leg parameters follow Eqs. (96) through (99), and 
by equating F 4 [Eq. (107)] to H B 1 4 , we find that 


NI 


B 4 


V* 



( 108 ) 
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As NI is reduced below NI Bi , ff 4 < H B , H 3 > H B , and H m remains above H g 
until NI reaches NI Bm , when H m becomes equal to H g . A^further reduction 
in NI to NI B3 lowers H 3 to H g , and if NI < NI g3 , then H < H g in every 
leg. Calculation of NI Bm and Nig 3 is more complex than calculation of 
N I B 4 , because of the dependence of some or all the leg parameters on the 
F values [Eqs. (94) and (95)]. 

To summarize, there are four regions of NI'. 

(1) NI > NI Bi , where H 3 > H . > H 4 > H B . 

(2) NI Bi > NI > NI b m , where N 3 > > H 4 . 

(3) NI Bm > NI > NI b3 , where H 3 > > H a > 

(4) NI b 3 > NI, where H g > H 3 > H m > H 4 . 

In calculating D versus the parameters of each leg will be assumed 

p and F 0 [Eqs. (88) and (89)] if H > H B , but if H < Hg, then the param- 
eters to be assumed are p' p and F' B [Eqs. (94) and (95)]. Computation of 
this type is extremely tedious when done manually, but is quite simple 
when performed on a computer. We notice, however, that in the region 
where NI > Nig 4 , all leg parameters are constant; hence, calculation of 
D in this region should be simpler, as we shall see next. 


d. D VERSUS (l 4 /lj) FOR NI > NI gi 

If NI > NI Bi [Eq. (108)], then the leg parameters p p and F 0 
[Eqs. (88) and (89)] are independent of NI . These leg parameters differ 
from each other only because of variation in dimensions, as expressed in 
Eqs. (96) through (99). Our objective now is to substitute Eqs. (96) 
through (99) into the previous derivations in order to express D versus 
leg-length ratios in a more direct fashion. 

The general solution, Eqs. (70) through (73), remains unchanged, 
except that the expressions for p and q are reduced to the following. 



(109) 
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and 



( 110 ) 


We shall now examine several special cases. 

When the core is unloaded, Eq. (110) is substituted into 
Eq. (74) directly in order to calculate D NL . 

For NI — <*> t substitution of the relation (p 4p /p 3 ^) * ( 1 3 / 1 ^ ) 
[Eq. (97)] into Eqs . (76) and (77) gives 


and 


PjVJ-co 


N 1 

h P _ *’ / *r _M 

*3 U, ' 0, / 

1 + — 


9/V/-0 


*3 

] 

<4 


I. 


1 + 


( 111 ) 


( 112 ) 


Substituting Eq . (72) into Eq. (Ill) results in a transcendental equa- 

tion in Pf/i-w- A numerical solution for Pp/j^a, can be obtained without 
any difficulty. The values of p NI ^ w and q NJ ^ [Eqs. (Ill) and (112)] 
are substituted into Eq. (73) in order to compute D N j _ a ,. For the no-load 

case, Eq. (78) becomes 


D 


NI , NI-™ 



,+ (i 3^4 


) ' 

7 . 


- S 


(113) 
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For the case of S 1, substitution of Eqs . (96) through (99) 

into Eqs. (79), (80), (81), and (81a) results in the following expressions 


»S-1 " 



(114) 


®NL , S- 1 


2 


2 - 



+ 


2 



(115) 


and 




1 ,NI^<x> 




D 


NL , S- 1 , AT/-® 


(116) 


(117) 


Note that only Eq. (117) yields the value of D that was expected when 
the problem of flux division was first introduced. 6 But one should 
recall that Eq. (117) was derived assuming no variation of p with flux 
in every leg, a condition valid only for small variations of 0 around 
0 a 0 in each leg, as was shown in Fig. 13. 

6. Summary 

An analysis of flux division ratio D in a loaded saturable core has 
been developed, assuming that w 3 = w 4 = ^ B /2 and that A0 r = 2 0 r . The 
calculation of A0 in a leg was based on a parabolic model for p. The 
final expression for D [Eq. (73)] is a function of three terms: 

5, which depends on the ratio of 0 r to 4 ) s \ which is a function of 
the ine 1 as t i c - swi tch ing parameters of each leg and SET MMF NI; and p, 
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which is a function of the leg parameters, /VI, and 4 ) r / < P Compu- 

tation of a general case involves a transcendental solution for p 
[Eqs. (70) through (72)]; if, in addition, NI is low, the six leg param- 
eters are determined by an iteration process that involves calculation 
of F in each leg and the nonlinearity of B p versus //. In the no-load 
case, the computation is simpler, because p = 0; hence, no transcendental 
solution for p is needed. The expression for D in the no-load case is 
given in Eq. (74), and if, in addition, NI -»cd, this expression is re- 
duced further, Eq, (78). It was shown in connection with variation of 
leg parameters with NI that, since A > A 7 = A , and Z 0 < l . , then 
3 > H m > // 4 during the whole switching time. The effect of l + / 1 3 on 
D is manifested indirectly by affecting the relative magnitudes of the 
leg parameters. Only in the no-load case, when NI is large enough to 
switch Leg 4 in its linear region, can the effect of l^/l^ on D be ex- 
pressed directly [Eqs. (74) and (110)]; e.g., Eq. (113) for D NL as NI — °°. 

The validity of this analysis will be checked in Sec. II-U by 
comparing computed and measured D versus NI for various load values. 

E. CONCLUSIONS 

The objectives in this section may be classified into two cate- 
gories: switching models (Parts A, B, and C) and the application of 
nonlinear switching models to the analysis of flux division (Part D) . 

In an attempt to improve our tools in dealing with switching prop- 
erties of a cons t an t - wi dth leg* the static <p(F) curve was calculated 
from assumed two hyperbolic B(H) functions, one for H £ h and the other 
for // >_ // th . These hyperbolic functions are superior to the common 
linear approximation of B(H) [ cf . Report 1, Fig. 4(b), p, 4], and were 
chosen simply from experimental observation. It is well recognized that 
analytic derivation of such models, based on some fundamental rules in 
physics, would be more elegant. However, until such models are derived 
analytically, having these empirical hyperbolic models is better for 
certain engineering purposes than having no models at all. 

Most of our switching models of a leg assume that the leg width, or 
rather cross-sectional area, is constant. Practical multipath cores, on 
the other hand, include variable-width legs. A typical example is the 
Transfluxor. Another example is the shaped MAD (see later. Fig. 42, p, 96) 
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are 


whose legs around the minor apertures, especially the inner ones 
variable in width. Calculation of switching parameters for a variable 
width leg of an arbitrary shape is very complex. For this reason we 
have simplified the problem by calculating p (a very important switch- 
ing parameter in common application) of a tapered leg as a function of 
geometry and flux level. Approximating half of the main leg in a 
Transfluxor, or Leg 3 in a MAD, by a tapered leg is more accurate than 
by a constant-width leg. Calculation of p of a tapered leg has, there 
fore, practical applications. 


Average p which is constant in time is useful when A 0 is calculated. 

In Report 1 (pp. 28-30), p and p ay of a constant-width leg were calculated 
only for full switching, i.e., forA0 = 2<t> r . In analyzing analog circuits, 
where <f> changes from ~<P r to an arbitrary level <p f , there is a need for 
the average p as a function of 0^.. We have shown that P ty ~ (1 /A0) P d/' 
and p = kcp/i^ 4, {\/p)< 14 >] and that p, not P ty , is the useful parameter for 
calculating A0. For this reason we have calculated families of p versus 
<P f /<P s , with <P r /4>, as a parameter for both a constant-width leg and a 
tapered leg. The point where p/p reaches maximum is the point where we 
can get the maximumA0 per unit of DRIVE excess charge- turns . 


No attempt has been made to use the expression for p in order to 
calculate flux division in a saturable core. However, inspection of the 
plot of ~p/p versus 4> f /4> s ( cf . Fig. 7) reveals the essential mechanism 

associated with flux division. Since l 3 ^ then Fq 3 < - Fq 4 and 

p i P > P 4 p : h ence > ^3 > ^<P 4 . This. i n turn, calls for P 3 /p 3p to be 

larger than P 4 /P 4p . which serves as another reason forA0 3 to exceed 
A (p^. Tt is this factor that complicates the dependence of D on l 4 /I 3 
and yields higher values of D than expected by simple considerations 
that ignore the parabolic relation between p and <p. 


We shall see in Sec. II-D how calculated D versus NT compare with 
measured data. This analysis will be extended in the future to a non- 
saturable core, such as a Laddie, once we have better switching models 
for a leg that is initially in a partially switched state. 
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II EXPERIMENTATION AND DISCUSSION 


The purpose of this section of the report is to verify experimentally 
and discuss some of the theory and calculation of Report 1 and Sec. I of 
this report. In Part A, some measurements of toroidal core switching 
properties are presented and discussed. Part B treats the design and 
geometrical considerations of two multipath cores that are to be used for 
experimentation on flux switching in multipath cores. The two final core 
designs, one saturable and the other nonsaturable, are also given in this 
section. In Part C, test procedures for obtaining the switching parameters 
of a leg are described in general, and an example is given using one of the 
designed cores. Experiments on flux division are described in Part D, 
using the two designed cores. The results for the nonsaturable core are 
compared with the published data on flux division in a Laddie (c/, Report 1 
pp. 75-80). The results for the saturable core are compared with the anal- 
ysis of flux division developed in Sec. I-D of this report. Finally, in 
Part E, the phenomenon of unsetting effect (c/. Report 1, pp. 96-99) is 
verified experimentally. 

A, PROPERTIES OF TOROIDAL CORES 

Most of the magnetic material properties can best be determined by 
choosing a simple geometry so as not to introduce superfluous complications. 
Thus, toroidal cores are used when material properties are to be studied. 

1. Effect of Switching Rate on Static 0(F) Curves 

One of the fundamental properties that is needed is the static 0(F) 
curve. The technique for obtaining this characteristic consists of using 
a positive dc MMF to set the core, and a negative rectangular pulse, super- 
imposed upon the dc MMF, to clear the core, as shown in Fig. 17(a) with 
Switch Sw open. The dc MMF acts during a time T between CLEAR pulses. 

Thus, for a low repetition rate (e.g., 60 cps), 0(F) is essentially the 

static 0(F) curve. The flux change is measured during the CLEAR pulse, 
since the high CLEAR MMF produces a 0 ( t ) of short duration and high ampli- 
tude, which is easy to integrate in order to determine A0. The setting 
begins at the time when the magnitude of the CLEAR pulse falls below the 
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(b) NET MMF APPLIED AND FLUX WAVEFORM 
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FIG. 17 MEASUREMENT OF STATIC 0(F) CURVE 
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dc MMF. The fall time of the CLEAR pulse is typically of the order of 
0.1 jJ.se c to 1 ji sec, so that . the dc SET MMF effectively comes on relatively 
fast. It was questioned whether this short rise time of the SET MMF 
significantly affects the resulting static 0(F) curve. There was no 
particular reason to expect a difference, but there was also no assurance 
that no difference existed. The following experiment was performed to 
see if there was any difference in the static 0(F) curve obtained with a 
short rise time compared with a very long rise time. 

The static 0(F) curve of a MAD, treated as a toroid, was determined 
for two different conditions. The first condition was the usual one de- 
scribed above. In the second condition, Switch Sw was closed, so that 
the 100-turn winding, which is in series with a diode, would constitute 
a short-circuited winding during the SET time, and an open -ci rcui ted 
winding during the CLEAR time. Alternatively, one could have increased 
the fall time of the CLEAR pulse, but this would have required a modifi- 
cation of the pulser. The current induced in this 100-turn winding slowed 
down the switching rate during the SET time by an appreciable amount. A 
comparison between the MMF’s and A0* s versus time in the two cases is 
shown in Fig. 17(b). An additional comparison of the two, cases is made 
in the 0-F plane in Fig. 18. The resulting two static 0(F) curves are 
shown in Fig. 19. These curves are in agreement with each other for 
values of F below the midpoint of the upper knee. Above this point, the 
curve for Sw closed is about 2.7% lower. It is not known presently 
whether this small difference is due to the difference in the switching 
rates, or due to an unaccounted-for experimental error. The extra wiggle 
in the steep part of Fig. 19 is a result of the existence of minor aper- 
tures in the MAD. It is concluded that the 0(F) curves taken with Sw 
open (without the diode) can be considered as being essentially static. 

2. Testing Blank Disks 

Special multipath cores, which were to be used in a study of flux 
division (to be described in Part D), were ul t rasoni cal ly cut from 54-inch 
diameter Telemeter Magnetics T-5 ferrite disks. To ensure that no exces- 
sive radial variations in material properties existed in these disks, 
they were individually tested as toroids of very high OD/ID ratio before 
the cutting process. This was done by recording a 10 fJsec 0(F) curve of 
each disk (each disk had a 0. 050- inch- di ameter hole in the center through 
which the SET and CLEAR windings could be wound). Disks having a sharp 
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FIG. 19 STATIC 0(F) CURVES WITH Sw OPEN AND CLOSED 

AMP Incorporated, Core No. 395813-1, Geometry is given in Fig. 42 









radial variation in coercivity, resulting in a 0(F ) curve with a kink in 
the steep part of the curve, were rejected. 

3. Modification of the Parabolic Model 

The parabolic switching model, 

4> = p(F - F 0 ) , (118) 

where 

P • p,[> -(^) 2 ] . < ii9 > 

agrees with experimental data at high F values, but is poor at low values 
of F [i . e . , F slightly greater than the lower knee of the static 0(F) 

curve]. This is evidenced both by the static 0(F) loop, which does not 

have vertical sides as is assumed in the model, and by the 0 p (F) curve, 
which is not linear at low F values, as is assumed by the model. It would 
be valuable to know for what reasons the model does not agree at low F 
values, at least phenomono logi c a 1 ly , and what approach one might use in 
modifying the model to obtain better agreement. For these reasons, ex- 
perimental 0(0) curves of a thin ring were examined foreach of several 
constant F values. Oscillograms of these experimental 0(0) curves can 
be obtained by the technique described in Appendix G. The thin ring had 
an OD/ID ratio of 1.06 and was made of a low- threshold ferrite (Telemeter 
Magnetics T-5). The 0(0) curve obtained from Eqs. (118) and (119) for 
constant F is parabolic and is, therefore, symmetrical in 0. The experi- 
mental 0(0) curves were symmetrical (except for the initial spike) at 
high F values, but very asymmetrical at low F values, with the peak to 
the left of the center of A0. One possible way of accounting for this 
asymmetry is to consider F 0 [Eq. (118)] to be some appropriate function 
of 0. F Q must increase as 0 varies from ”0 r to 0 4 in order to shift the 
peak of the 0(0) curve to the left. Such an F Q (0) variation also results 
in a static 0(F) loop whose sides are not vertical, but have finite posi- 
tive slopes. This is observed experimentally even for very thin rings. It 
is possible that the effects of radial variations of H resulting from an 
OD/ID ratio significantly greater than unity can also be accounted for by 
an appropriate F Q (0) function. Since the final value of 0 attained for 
each constant value of F is described by the static 0(F) curve, it is 
likely that the inverse of the static 0(F) curve is a reasonable approxi- 
mation for F q ( 0), especially toward the end of the switching. 
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An attempt was made to account for the Fg(0) variation by assuming 
that F Q is equal to d of the static 0(F) curve (c/. Report 1, Fig. 9, 
p. 8). It is also assumed that the initial value of 0, 0 Q , was equal 
to ~0 r , The resulting calculated 0(0) curves, which contained a constant 
P p , yielded a nonlinear $> (F) function that was still different from the 
experimental 0^(F) curve. Therefore, to obtain a reasonable comparison 
of computed and experimental 0(0) curves, p p was assumed to be nonconstant 
with F and its value calculated for each value of F by equating the com- 
puted and experimental 0^ . The resulting calculated and experimental 0(0) 
curves are shown in Fig. 20 for three values of F . The values of p p are 
also given in Fig. 20. The 0(0) curves are plotted normalized to 0 f , whose 
value is obtained from the experimental static 0(F) curve [Fig. 21(a)] at 
5F q , or F 5 7.5 ampere- turns. Also shown in Fig. 21(a) are the values of 
F corresponding to the 0(0) curves of Fig. 20. This static 0(F) curve is 
the one that was used for the F Q (0) variation in computing the 0(0) curves 
of Fig. 20. Note that if F Q had been assumed constant [i.e M the static 
0(F) loop had vertical sides] , then the 0(0) curves of Fig. 20 would all 
terminate at 0/0 f B 1. The experimental 0^(F) curve for this core is 
shown in Fig. 21(b). Note that the values of F in Fig. 20 are in the very 
nonlinear region of this <p p {F) curve. The initial spike observed in the 
calculated 0(0) curves of Fig. 20 is a result of the low values of F Q ob- 
tained from the static 0(F) curve near 0 - -0 [see Fig. 21(a)], 

It is concluded that the parabolic model can be significantly im- 
proved for low F values by considering variation of p p with F and varia- 
tions of F 0 with 0, but that some disagreement still persists. The next 
step in the modification is the determination of analytical functions for 
P p (F) and F o (0). The major disadvantage of the modification of the para- 
bolic model is, of course, the added complexity, which in some cases may 
not be justified. It may be helpful for some cases to introduce P p (F ) , 
but not F Q (0) . In this case, p p (F) can be determined from the0^(F) curve, 
as was done in Sec. I-C. 

4. Summary 

The effect of switching rate during the setting time of the static 
0(F) experiment was found to be very small ( cf . Fig. 19). Blank disks, 
from which multipath cores were cut (for use in the f 1 ux -d i vi s i on experi- 
ments), were selected on the basis of their 0(F) property. Experimental 
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FIG. 20 COMPARISON OF EXPERIMENTAL AND COMPUTED #/>) CURVES 
FOR FERRITE THIN-RING E-6 

r ( * 3.53mm; r Q * 3.75mm; r av = 3.64mm; r Q /r. = 1.06; h ■ 0.688mm 
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(o) STATIC </>(F) CURVE 



FIG. 21 OPERATION IN FIG. 20 RELATIVE TO STATIC AND DYNAMIC 
CHARACTERISTICS 
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0(0) curves of thin ferrite ring are compared with curves calculated from 
the parabolic model (Fig. 20), with the following modifications: The F Q 

of the parabolic model is assumed to be a function of 0. This F^ (0) func- 
tion is approximated by the inverse of the static 0(F) curve. The values 
of p in the parabolic model were assumed to be some function of F and are 
determined by equating calculated and measured peak 0 values, 

B. GEOMETRY CONSIDERATIONS OF MULTIPATH CORES 

As discussed in Report 1, p. 1, a multipath core can be divided into 
legs; each leg is characterized by appropriate switching models. Assuming 
uniform material, the switching properties of the legs differ from each 
other only due to variations in geometry. Unfortunately, the geometry of 
a leg is sometimes not well defined. We wish to discuss this problem in 
connection with the design of two multipath cores that were used for ex- 
perimental verification of phenomena associated with core geometry in 
general, and flux division in particular. The design and analysis of these 
two cores involved a number of geometrical considerations that will now be 
discussed . 

1. Design Considerations 

а. SATURABLE AND NONSATURABLE CORES 

Two three-leg cores were designed for the f lux- d i vi s ion ex- 
periment. The first core was intended for use in measuring flux divi- 
sion in a nonsaturable core. With this core, we checked and supplemented 
the experimental data on another core (Laddie) that was reported in the 
past. 6,7 We also hope to use these data for comparison with future 
analysis. The second core was saturable and served to verify experi- 
mentally the analysis of flux division developed in Sec. I-D. Both cores 
were u 1 tras on i c a 1 ly cut from ferrite disks of constant thickness. 

б. CROSS-SECTIONAL AREA OF A LEG 

One of the most important geometrical properties of a leg to 
consider is the variation of its cross-sectional area or, if the thick- 
ness is constant, the variation of its width. There are many problems 
associated with the description of flux switching in a va ri ab le- wid t h 
leg. These problems were discussed in Sec. I-B-2. Because of these 
problems, the experimental work on flux division was performed on 
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multipath cores that were spec i a 1 ly designed so that each leg had constant 
width for as much of its length as possible. Once flux division for such 
cores is well understood and predictable, there is a possibility of under- 
standing flux division in less ideal cores. 

C . JUNCTION PROBLEMS 

Ideally, the legs of a multipath core should have dimensions 
that are uniquely defined, irrespective of the switching paths in the 
core. These switching paths are determined not only by the location of 
the drive windings, but also by the amplitudes of the drive currents. 

Each leg should have a well defined short edge of length Z^, ( cf . Fig. 19, 
Report 1, p. 18), long edge of length Z qI and cross section of area A . 

In addition, the midpaths of all legs should meet at nodes, or junction 
points, from which the values of Z aT may be obtained [actually for very 
thick legs, a path somewhat inside of the midpath should be chosen, 10 
because of the difference between l/(Z aT ) and (1/I)* T ]. These ideal con- 
ditions are assumed by the leg model, but cannot be realized physically. 
However, by designing cores that minimize the junction problems and by 
making suitable approximations in the determination of leg dimensions, 
the leg model can still be very useful. The alternative to using the leg 
model, i.e., solving the non-linear field problem, is a very complex 
problem. In order to minimize junction problems, the junction volume 
should be made as small as possible, so that the percentage of a leg 
dimension that is ambiguous is minimized. This is done primarily by 
making the ratio of leg length to width, l/w , as high as possible. 
However, a high l/w ratio has the problem of flux leakage when a lumped 
drive winding is used. As a compromise, a value of l/w ^ 18 was chosen 
for the longest leg. Fillets should be avoided, since they also increase 
the junction volume. In designing a multipath core for which all legs 
have the same width, there are two advantages in having the three legs 
radiate from the junctions at approximately equal angles from each other: 
First, flux leakage is minimized; second, the location of the junction 
node is less dependent on where switching takes place. 

d. LEG CONTOUR 

The presence of sharp inside corners in a leg is not desirable 
because of leakage flux that bypasses the corner (Report 1, Fig. 106, 
p. 109). For this reason, unnecessary sharp corners were eliminated in 
the core designs. 
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e . FINAL CORE DESIGNS 

The nonsaturable core, Core N, and the saturable core, Core 5, 
were designed on the basis of the above discussion. The dimensions of 
these cores are shown in Fig. 22. Note that Core N has three legs of 
the same width, whereas Core 5 has two legs of the same width and a third 
leg of twice that width. 

2. Determination of l (t l Ql and l my 

In order to' obtain the switching parameters of the legs in the cores 
under discussion, we must first determine the dimensions l Q , and 

of each leg. However, the exact dimensions of each leg depend upon the 
paths of switching. To illustrate the problem, consider Core N and com- 
pare the case when the drive is on Leg 1 with the case when the drive is 
on Leg 2. The values of I. and l Q of each leg are different in these two 
cases, as shown in Fig. 23. These particular choices of dimensions in- 
volve compromises and approximations, some of which will be discussed in 
the next part. 

3. Flux Switching in a Core Having a Re-entrant Shape 

Discussion of the topic of this section requires the use of termi- 
nology that first must be defined. 

A switching line in a magnetic core is a line at every point of 
which the vector M is tangential. A switching path in a magnetic core 
is a group of switching lines. 

The leg length l " 1 n is determined by the line along which inelastic 
switching first occurs in a cleared leg (<fc *= ~4> r ) as the MMF is slowly 
increased from zero. The leg length Z max is determined by the line along 
which inelastic switching last occurs as the MMF continues to increase 
slowly. Thus, [® in and l ro a x are determined by lines that may or may not 
follow the short and the long edges of the leg, as will be illustrated next. 

Consider a core having a re-entrant shape, such as obtained by removing 
Leg 1 of Core N [cf. Figs. 22(a) and 24]. Assume that this core is ini- 
tially well cleared in a clockwise sense, and that a dc current is applied 
to the winding on Leg 3 as shown in Fig. 24. The magnitude of this current 
is slowly increased until switching just begins in Leg 2. The initial 
switching line will be the shortest closed line linking the drive winding, 
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FIG. 22 DESIGNS OF CORES FOR EXPERIMENTATION 

Dimensions in mm; h = 1.31mm; Tolerance * ±0. 010mm 
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FIG. 23 DETERMINATION OF /. AND l o (in mm) FOR 
EACH LEG OF CORE N 
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Consider Point E in Fig. 24. Since the core was initially well 
cleared, the magnetization M followed the contour of the core. Thus, the 
initial magnetization at Point E can be represented by the vector M- in 
Fig. 24. When the MMF is increased just enough to establish switching in 
Line A-B-C-D-E-A, M changes from initial value M i to an intermediate 
value The vector differenceAM is parallel to M, which is tangential 

to the switching Line A -B -C-D -E -A . Whenever M at a point in a magnetic 
material changes direction by an angle other than 180°, some transverse 
switching is occurring and M at that point is non-parallel to M. In the 
absence of bulk anisotropy, this transverse switching must be accompanied 
by a transverse H field. The source of such an H field is, in general, 
made up of two components: the applied H, and the H due to div M (effec- 
tive magnetic pole density), which develops on the surfaces and inside the 
magnetic material (see Report 1, p. 128; also Ref. 11), The properties 
of transverse switching need further investigation. 

The above description is intended only to present the problems in 
understanding the switching in a core, such as the one of Fig. 24, and 
to estimate the location of the switching lines for slow switching. If, 
instead of the slow switching, fast switching is produced as a result of 
applying a relatively h i gh - amp 1 i t ude rectangular current pulse, then dif- 
ferent switching lines would be obtained because switching would be oc- 
curring simultaneously, though at different rates, throughout the entire 
core. Several other examples in which the re-entrant shape is also en- 
countered are the following: 

(1) Driving on either Leg 2 [see Fig. 23(b)] or Leg 3 of Core N 
involves switching in the series combination of Legs 2 and 3, 
which has a re-entrant shape. 

(2) In determining Z i3 inside the junction of Fig. 23(a), the 
problem is mixed with the junction problem. 

(3) The re-entrant shape also occurs in Legs 2 and 3 of the multi- 

aperture core shown in Fig. 97(a), (p. 102) of Report 1. 

The detailed process of switching in a core as in Fig. 24 is diffi- 
cult to determine, even for slow switching. However, if the shape of the 
core is made more re-entrant, as in the core shown in Fig. 25, and the 
material is assumed to have a static B(H) loop with vertical sides, then 
a simple model for slow switching can be proposed. This model has the 
zone pattern shown in Fig. 25, which provides magnetization closure and 
crossover switching paths ( i.e. f switching paths which cross over from 
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FIG. 25 SWITCHING PATH CROSSOVER ZONES IN 
A CORE HAVING A RE-ENTRANT SHAPE 
Arrows Indicate Directions of M 

one edge of the leg to the other edge). For a closer examination of these 
zones, refer to Figs. 26(a), (b) and (c). The dashed lines in these fig- 
ures represent the zone boundaries. Note that the relative directions of 
the magnetization M and the zone boundaries are such that no magnetic poles 
exist at any of the zone boundaries (c/. Report 1, p. 122). The switching 
paths are identified by arrows in Fig. 26(c). Each switching path is bound 
by a dashed line and a dash-dot line. Note that the switching path could 
not pass diagonally through the rectangular zone without establishing new 
zone boundaries within the existing zone. Figure 26(c) also illustrates 
the motion of the zone boundaries. Consider the three points, Pj , P 2 , and 
P 3 in Fig. 26(c). The initial M, , the final M, , and several inter- 

mediate values of M at each of these points are shown in Fig. 26(d). The 
physical mechanism of the switching required at Points P 2 and P 3 has not 
been thoroughly investigated. The vector at Point P 2 is parallel to 
M . at Point P^. Note that as the zone boundary moves past Point P^, M 
switches by 180°, whereas at Points P 2 and P 3 , M switches by less than 
180°. Material at a point near the rectangular zone switches first as 
at Point P 2 and later finishes switching as at Point P 3 , resulting in 
a net 180® switching after both zone boundaries have passed the point. 





(c) SWITCHING PATHS 


M ( P, M, 

— • ► 


Am 


(d) SWITCHING OF M AT SEVERAL POINTS IN (c) 

RB- 3696-1 55 

FIG. 26 SWITCHING PATH CROSSOVER ZONE 
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Also notice that the vectors AM, which are parallel to M have the same 
directions as the zone boundaries. Based on the definition of a switching 
path, these zone boundaries coincide with the switching paths. The zone 
boundaries at a later time in the switching process are compared in Fig. 27 
to the ones obtained earlier. The angle 6 that the boundary of the rectan- 
gular zone makes with the edge of the core will be determined for a partic- 
ular core shape by minimization of the total magnetic energy of the system. 
This includes primarily z one - bound a ry energy and magnetostatic energy, 
which results from both the applied H and the H due to the existence of 
magnetic poles. For different angles 8, the relative widths of the 
switching paths near Points P 2 and P 3 in Fig. 26(c) will be different. Also, 
the angle between the initial and final direction of M at Points P 2 and P 3 
in Fig. 26(d) will be different. These two effects change with 8 in such 
a way as to provide the required continuity of 4> in the switching path, 
and also to divide <£> equally between the two branches of the switching 
pat h . 
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FIG. 27 COMPARISON OF CROSSOVER ZONES EARLY 
AND LATE IN SWITCHING 


The extension of this zone model to materials whose static B(H) loops 
do not have vertical sides results in a configuration in which all of the 
material between two zone patterns (the shaded regions in Fig. 27) is in the 
process of switching. These two zone patterns then represent the initia- 
tion and completion of switching. 

This switching path crossover zone model has not been verified ex- 
perimentally. It represents a first attempt at describing the details 
of crossover switching. 
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4. Machining of Multipath Cores 

The multipath cores were u 1 1 r ason ica 1 ly cut from 0.050-inch thick 
ferrite disks. The ultrasonic tool was machined with a three dimensional 
pantograph from an enlarged (32*) plastic and metal pattern. This tech- 
nique, although quite workable, has several problems. The major problem 
is the development of a taper in the cross section of a core leg. This 
problem is not easily remedied, and is most pronounced for very thin legs. 
Another problem stems from inaccurate dimensions. A leg width can easily 
be off by 0.05 mm, unless special care is taken in designing and construc- 
ting the ultrasonic tool. Cores N and S, used in the experiments of this 
report, have dimensions that are within ±0.010 mm of the design dimensions 
shown in Figs. 22 and 23. This accuracy was accomplished with consider- 
able difficulty. 

5. Summary 

The various considerations involved in designing two cores for ex- 
periments of flux switching in multipath cores are discussed. The most 
significant of these considerations is the junction problem. The two 
final core designs are shown in Fig. 22. 

The problems associated with determining the leg parameters 1 ; , l 0 , 

and i* T are considered, and an illustration of the different results 
corresponding to two different drive conditions on Core N (Fig. 23) is 
given . 

The terms switching line and switching path are defined, and then 
used in examination of switching on cores having re-entrant shapes. 

Only a crude picture is given for switching in one core (Fig. 24), where- 
as a model is proposed for describing switching in a second core (Fig. 25). 

C. DETERMINATION OF LEG PARAMETERS 

A number of leg parameters must be determined for each leg of a 
multipath core so that switching models can be used to describe the 
switching phenomena (such as flux division). If the material parameters 
are known, these leg parameters can generally be calculated from the leg 
dimensions; however, typically there are significant variations in mate- 
rial parameters from core to core (and even from point to point within 
a core). Therefore, it is desirable that tne leg parameters be determined on the 
same core that is used in the experiments, 'lhe particular disks from which 
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Cores JV and 5 were cut are believed to have had some radial variations 
in threshold* It is difficult, however, to measure leg parameters sepa- 
rately for each leg of a core. As a compromise, the material parameters 
were experimentally determined for two legs in series (comprising the 
closed path needed for the experimental measurements); these same material 
parameters were then used to calculate the leg parameters for each leg of 
of the same core. Material variations from core to core were thus ac- 
counted for, but material variations within a core were neglected. 

This section describes experimental techniques for determination of 
the required leg parameters. The details of the experimental measurements 
are included in Appendix H. 

1. Parameters to be Determined 

A list of leg parameters is given in Table I, together with the mate- 
rial parameters from which these leg parameters are determined and the ex- 
periments from which the material parameters are determined. 


Table I 


DETERMINATION OF SWITCHING PARAMETERS 


EXPERIMENTS 

MATERIAL 

PARAMETERS 

LEG PARAMETERS 
OR PROPERTY 
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flux switching 
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4> (F) 

p 
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2* Closed Path for Testing 

The basic experiments listed above must be performed on a closed 
magnetic path. Thus, some closed path in the multipath core must be 
chosen and flux switching in any other paths must be blocked. One way 
in which this blocking can be accomplished is by putting a short-circuited 
winding around the legs in which no switching is wanted. For example, the 
parameters of Core N [Fig. 22(a)] were determined by choosing Legs 1 and 2 
for the closed path and putting a short-circuited winding on Leg 3. 
Similarly, for Core S [Fig. 22(b)], Legs 3 and 4 were chosen for the closed 
path and Leg m was linked by a short-circuited winding. If possible, this 
closed path should have the same cross-sectional area over its entire length 
otherwise the determination of material parameters would become more in- 
volved and less accurate (c/. Report 1, pp. 109-112). 


Blocking 4> in the legs where no switching is wanted may present a 
problem because of the resistance and leakage inductance of the short- 
circuited winding. The effectiveness of a short-circuited winding for 
preventing switching in a leg, if the leakage inductance of the winding 
is neglected, is a function of N 2 /R, where N is the number of turns of 
the winding and R is its internal resistance. Assuming that the avail- 
able aperture is filled with the short-circuited winding, R is approxi 
mately proportional to N 2 , so that the effectiveness of the short-circuit 
is not significantly dependent upon the number of turns used in the winding. 
On the other hand, if the resistance of the winding is assumed to be zero 
and the effect of leakage inductance of the winding is calculated, the 
pertinent quantity is N 2 /L. The leakage inductance is approximately of 
form N 2 K,[ln(K 2 SN)~K 3 ], where K , , K , . and K , are constants, so that a 
winding of one turn that fills the aperture is better than a multiturn 
winding that also fills the aperture. Additional factors in favor of a 
single turn are (a) a better filling factor, and (b) the possibility of 
using a copper strap, heavier than the wire inside the aperture, in order 
to complete the short-circuited turn and thereby lower the resistance and 
leakage inductance. The short-circuited turn might also be made by immer- 
sing the leg to be blocked in mercury or some other liquid conductor. 

When performing the experiments, it is a good idea to check for flux 
switching in the blocked leg by putting a sense winding on it. Another 
possible means of blocking a leg in which no switching is wanted is the 
physical removal of the leg. This means requires either that the flux 
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division experiments be done before the leg is removed for the determi- 
nation of leg parameters, or that two separate cores are used, in which 
case they must be known to have the same magnetic properties. 

3. Static 0(F) Experiment 


The static 0(F) experiment consists of two measurements: 

(1) Measurement of the total flux change (inelastic plus 
elastic) from negative remanence to any intermediate 
level up to positive saturation. 

(2) Measurement of the elastic flux switched from negative 
remanence toward negative saturation. 

In the first measurement, the core is set with a pos i t i ve dc MMF and cleared 
with a negative MMF pulse whose amplitude is sufficient to attain a repro- 
ducible residual state. The flux change is measured during the CLEAR pulse, 

# 

since the high CLEAR MMF produces a 0(t) of short duration and high amplitude, 
which is easy to integrate. In the second measurement, the dc MMF is removed 
and the amplitude of the CLEAR pulse is varied. The test procedure for ob- 
taining the static 0(F) curve is outlined in Appendix H. 


The resulting static 0(F) curve is then used to determine the material 
parameters B r , B ss , H a , H n , and H ( cf . Fig. 2, Sec. I-A). First, B r is de- 
termined from 


<K_ 

hw 


( 120 ) 


Next, Eq. (16) is used to determine B ss and // from the second experiment. 
This. can be done either by using two experimental points, or by the method 
of least mean squares using all of the experimental data. For most appli- 
cations, Eqs . (17) and (18) can be used instead of the more complicated 

Eqs. (A-6) and (A-7) for the determination of H and H . This can be done 

n <7 

either by using two appropriate experimental points, one for Eq. (17) and 
the other for Eq. (18), or by the method of least mean squares using all 
of the experimental data from the first static 0(F) curve. 

Once the material parameters B ss , B r , H a , H n , and H are determined, 
"th is calculated from Eqs. (12) and (13), and then the static 0(F) curve 
can be computed for any leg of the same material by applying these param- 
eters and the leg dimensions to Eqs. (16), (17), and (18). 
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As an example, the computed B(H) curves of Core N are shown in 
Figs. 28(a) and (h), and expe r imen t a 1 and computed static 0(F) curves of 
Legs 1 and 2 of Core N in series are compared in Figs. 28(c) and (d). 

The values of the material parameters are given in Fig. 28. 

The agreement between the measured and the calculated static 0(F) 
curves is very good, except where 0 < F < an ^ where \f\ > 10 

The agreement for 0 < F < 0 . 37 AT i s not good because of some inelastic 
flux switching that is neglected by the model. For |F| > 5 AT, the com- 
puted curve (even if 0 # . r is included) has poorer agreement with the ex- 
perimental data than for |f| ^ 5 AT, If the experimental data were fitted 
for the range of values 100 AT < F < 0 , then the value of & ss obtained 
would be larger than the value required to obtain good agreement for 
.1, < F < 10 H . I .. We thus conclude that for the range of F values 
most frequently used in computer applications, the static 0(F) curve may 
well be described by Eqs. (16), (17) and (18). 

The material parameter B $ is obtained from the experimental static 
0(F) curve by choosing the value of 0 5 to be at about five times F^. Thus, 


B 

* 




(121) 


From this value of B s , one can calculate the <p , value for any leg by using 
the h and w of the leg. 

4. 0(F) Experiment 

The 0 p (F) experiment consists of setting and clearing the core with 
rectangular pulses and measuring the peak 0 as a function of the SET pulse 
amplitude. The test procedure is outlined in Appendix H. As discussed in 
Sec. I-D-5, the resulting experimental 4> p (F) curve is described by Eq. (87) 
in the high-F linear portion, and by Eq . (90) in the low-F nonlinear 

portion (c/. Fig. 15). The nonlinear approximation defines k, F", and v. 
Which of these approximations is used depends upon the accuracy needed, 
the range of F for which the parabolic switching model is to be used, and 
the degree of complexity one can handle for a particular application. 

From these experimental parameters (which apply to the entire closed mag- 
netic path) the corresponding material parameters can be calculated by 
using Eqs, (88) and (89) or Eqs. (91) and (92). These material parameters 
can then be used to calculate the leg parameters for any other leg of the core. 
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FIG. 28 CALCULATED AND MEASURED STATIC B(H) AND 0(F) CURVES FOR 
LEGS 1 AND 2 OF CORE N IN SERIES 

B s , = 0.2270 wb/m 2 ; B, - 0.21 14 wb/m?; 0 r = 0.2567 • 10"6 wb; 

H n = 32.0 A/m; H q « H, h = 36.0 A/m; H a = 230 A/m; 
l\\2 ~ 10.56mm; / 0 23 = 16.47mm; w = 0.927mm; h = 1.31mm. 
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(d) <f> (F) IN NEGATIVE SATURATION 


FIG. 28 Concluded 
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A numerical example can be found in Fig. 35, where measured and calcu- 
lated 4>?(F) of Legs 3 and 4 of Core 5 are compared. 

5. Summary 

The parameters which need to be determined for each leg were sum- 
marized in Table I, together with the corresponding material parameters 
and also the two experiments used to determine these material parameters 
The requirement of choosing a closed path in the core to be tested and 
the problem of preventing flux switching in all other paths were treated 
Both of the two experiments, the static 4>{F) experiment, and the <fi p (F) 
experiment, are presented and the process of using the experimental data 
to obtain the parameters for any leg of the same core is discussed. 

D. FLUX DIVISION 

1. Introduction 

In dealing with flux division, we distinguish between a saturable 
and a nonsaturable multipath core. As explained in Sec. I-D-l, all the 
legs of a saturable multipath core can be brought simultaneously into 
saturation residual state of <p = ±0 r , whereas in a nonsaturable core, 
at least one leg is unsaturated. Flux division in a saturable core was 
analyzed in Sec. I-D. Due to the lack of switching models for a leg 
initially in a pa r t i a 1 1 y -swi t ched state, flux division in a nonsaturable 
core has not been analyzed. Experimental data on flux division in a non- 
saturable core (the Laddie) was described when the problem of flux divi- 
sion was first introduced. 6 

In this part, we first show results from extending the experimental 
investigation of flux division in a nonsaturable core. Second, verifi- 
cation for the analysis of flux division in a saturable core (Sec. I-D) 
is given by comparing measured and computed plots of D versus NI for 
various load levels. Along with this verification, the effects of 
various factors (such as 0 r /0 jf leg parameters, load levels, and Z 4 / Z 3 ) 
on D are investigated. 

2. Flux Division in a Nonsaturable Core 

a. experimentation 

The nonsaturable core used in the experiments on flux division 
was Core N [Fig. 22(a)], whose three legs are equal in cross-sectional 
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area, but different in length. As explained in Part B-l, this three-leg 
core was especially designed and fabricated for the experiment on flux 
division. The manner of dividing Core N into legs and the lengths of 
these legs are given in Fig, 23. 

In order to obtain a reproducible initial state, Core N was 
heated to slightly above its Curie temperature for about two minutes. In 
this way, the entire core was brought into a uniformly demagnetized state 
before any measurement was performed. 

First, MMF NI was applied to Leg 1 [Fig. 23(a)], and the re- 
sulting c£j , 0 2 , and were sensed and integrated. With no load on either 
leg, the experimental setup was similar to the one shown in Report 1, 

Fig. 72(a), p. 75. The bipolar SET and CLEAR current pulses were equal in 
amplitude and duration (5 microseconds). As NI was mono t oni c a 1 ly increased 
from zero, residua 1 A0j , A0 2 , andA0 3 were measured as a function of NI , 
using the f lux- re f e rence technique (c/. Appendix F). The same procedure 
was then repeated with the NI drive applied to Leg 3. 

6. RESULTS 

Plots of measured A0J , A0 2 , and A0 3 versus NI are shown in 
Fig. 29. Figure 29(a) is for NI applied to Leg 1, and Fig. 29(b) is for 
NI applied to Leg 3. Corresponding plots of flux division ratio versus 
NI are shown in Fig. 30. 

Oscillograms of <f> i , <t> 2% 4> z and NI for the case of drive on Leg 1 

are shown in Fig. 31 for NI - 0.7 AT, and in Fig. 32 for NI - 3.0 AT, 

These oscillograms were taken at a temperature of 29°C and repetition rate 
of 500 cps. The waveforms of 0j, 0 2 , and 0 3 during the SET time for the 
case of NI 55 3.0 AT are shown in Fig. 32(a), an enlargement of during 

the SET time is shown in Fig. 32(b), and the waveforms of 0j , 0 2 » an( * ^3 

during the fall time are shown in Fig. 32(c). 

These experimental data will now be discussed. 

C. DISCUSSION 

(1) In examining Fig. 30, the inaccuracy of theA0 measurement 
should be considered. The flux change in each leg was 
measured by means of the flux reference (c/. Appendix F), 
whose error is less than ±0.1 maxwell turns. The lower NI 
is, the smaller is the measured Acp in the longer leg; hence, the 
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(b) DRIVE ON LEG 3 


FIG. 29 MEASURED FLUX DIVISION IN A NONSATURABLE CORE 
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FIG. 30 FLUX-DIVISION RATIO IN NONSATURABLE CORE N 



TIME — > 


R P— 3 6 9 6— 197 

FIG. 31 FLUX-DIVISION OSCILLOGRAMS OF CORE N FOR 
Nl = 0.7 AT ON LEG 1 

Time scale = 1 /xsec/major div.; 
c p scale = 0.09 (V/T)/major div. 
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(a) SET TIME 

Time scale = 0.1 yusec/major div. 
4 > scale = 1.6 (v/T)/major div. 




R P- 3696-196 

FIG. 32 FLUX-DIVISION OSCILLOGRAMS OF CORE N FOR Nl = 3.0 AT ON LEG 1 
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higher is the error. Thus, when NI is applied to Leg 1, 
at NI = 0.592 AT, A0 2 /A <£ 3 reaches a peak of 82.3 ^-\\\\ 

(i.e., the error is less than - 20 ’.& however, for 

NI = 4.138 AT, (A0 2 /A0 3 ) = 11.91 ±0.16, (i.e., the error 
is less than ±1.34%). When NI is applied to Leg 3, at 
NI = 0.794 47\ A<£ 2 /A<£j reaches a peak of 10.7 I 
(i.e., the error is less than -^'7 %) I however, for 
NI = 2.914 AT, (A0 2 /A0j) = 3.08 ± 0.12, (i.e., the error 
is less than ±3.9%). These errors for the low values of 
NI could have been decreased by increasing the number of 
turns of the sense windings on the longer legs. This im- 
provement was made in measuring the flux division in Core S, 
as will be shown later. 

(2) It is hard to compare our experimental data with the data 
in Refs. 6 and 7 because we do not know what material was 
used for fabrication of the Laddie reported on in Refs. 6 
and 7. There is, however, similarity between the plots of 
Ac£ 2 andA <£ 3 for the ratio of l 3 /l 2 = 3, which is common to 
Refs. 6 and 7 and Fig. 29(a), except for the following. 

In Refs. 6 and 7 [see also Report 1, Fig. 72(b), p . 7 5 ] , A <fi 2 
and Ac /> 3 versus NI is less curved and rises with a higher 
slope than in Fig. 29(a). This difference might have been 
caused by the use of material of higher threshold and 
squarer B(H) loop in Refs. 6 and 7. Another difference be- 
tween these sets of data is that in Fig. 29(a), A0 2 and A0 3 
versus NI have small but finite slope at high values of 
NI (say, above 2 AT), whereas in Refs. 6 and 7, A <^ 2 and A ^ 3 
are essentially constant in this range of N I . 

(3) In the plots of D versus NI for leg ratios of 3 and 1.5 
(Fig. 30), we observe that in each core, D reaches a peak 
at some relatively low value of N I . In order to understand 
what causes this peak to occur, let us refer also to 

Fig. 29(a), where NI is applied to Leg 1 . As NI is increased 
from zero, there are slight inelastic flux changes in all legs 
until at around NI e 0.4 AT, H in Legs 1 and 2 reaches the 
threshold value. A further increase in NT causes an appre- 
ciable inelastic A0 in Legs 1 and 2, but A0 3 hardly increases, 
since // 3 is below the (soft) threshold. In this region of 


81 



NI , there fore , A 0 2 /A<£ 3 increases with AT. As AT exceeds 
''-0.6 AT, II 3 exceeds the (soft) threshold, and A0 3 increases. 
For values of AT above 0.8 AT, A0, andA0 2 start to level 
off due to saturation, but A<£ 3 continues to rise. This 
calls for a decrease inA 0 2 /A0 3 . Thus, at a value of AT 
somewhere between 0.6 AT and 0.8 AT (Fig. 30), A0 2 /A0 
reaches a peak. A similar mechanism causes the peak in 
A0 2 /A0, versus AT in the case where NI is applied to Leg 3 
except that since l 3 > l^ t this peak is reached at a value 

of NI higher than when NI is applied to Leg 1. 

(4) The oscillograms of 0 2 and 0 3 , (Figs. 31 and 32), are 

similar to the ones in Ref. 7, except that in Ref. 7 these 

oscillograms were not extended to the fall time. As in 
Ref. 7, 0 3 becomes negative (while AT is still on) if NI 
is high enough. The explanation for this phenomenon can 
be found in Report 1, pp. 77-80. 

(5) The waveshapes of 0, , and ^ during the fall time of 

AT were predicted in Report 1 [Fig. 75(b), p. 79]. Com- 
paring peak values of 4> in Report 1 , ji 3p I > \d> 2p \, whereas 
in the oscillograms of Fig. 32(c), |0 3p | < J 0 |. This 

disagreement stems from the following reason: Prediction 
of the waveforms was based on the assumption that since 
Leg 3 is in a soft state and Leg 2 is in a hard state, the 
slope e - (<14>) e /dF in Leg 3 is higher than in Leg 2. Such 
an argument is correct but incomplete; there is an addi- 
tional factor, i . e . , Leg 3 being longer than Leg 2, which 
acts in an opposite direction. Which factor is predominant, 
depends on the circumstances. Apparently, the factor of 
leg-length ratio is predominant in Fig. 32(c), whereas in 
Fig. 75(b) of Report 1, the soft-state factor was assumed 
predominant. Based on this point, and since (dd>) e /dF 
varies during the fall time, instead of comparing^,, and 
and <P p3 , we should have compared A0 f2 andA0 f3 during the 
fall time. When this comparison is made in Fig. 32(c), we 
find that |A0 e3 | > |A<£ e2 |, which agrees with the basic as- 
sumptions made in Report 1 in predicting the waveshapes of 
in the three legs. 
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3. Flux Division in a Saturable Core 

а. EXPERIMENTATION 

Core S [Fig. 22(b)] was designed and fabricated specially for 
measurement of flux division in a saturable core. The manner of dividing 
this saturable core into Legs m, 3, and 4 and the dimensions of these legs 
are also shown in Fig. 22(b). 

The experimental setup followed the one shown in Fig. 10, p. 26. 
Leg 4 was first unloaded. The core was cleared to a negative remanent 
state (not shown in Fig. 10). A current pulse of amplitude I and duration 
T was then applied to an /V-turn winding on Leg m in the SET direction, and 
the net flux change in Leg m, ^4 ) m , was measured. The magnitudes of I and T 
were adjusted until was equal to 2 4> ft which is the maximum possible 

inelastic A0 in Leg 3 (or Leg 4). For each pair of values of I and T that 
yielded A^ ■ 2<£ r , the values ofA$ 3 ,Ac£ 4 , and I were recorded. 

This procedure was repeated for three values of load: 

(1) N l = 2, R l = 1.02 ohm, i.e., (,N 2 L /R L ) = 3.926 T 2 /ohm. 

(2) N l = 1, R l = 0.093 ohm, i.e., (N 2 L /R L ) = 10.66 7’ 2 /ohrr. 

(3) N l = 2, R l = 0.1 ohm, i.e., (N 2 L /R L ) = 40.00 r 2 /ohm. 

б. RESULTS 

The resulting plots of A<£ 3 and A<£ 4 versus Nl are shown in 
Fig. 33 for the four values of N*/R L . The corresponding plots of D will 
later be compared with calculated plots. 

Oscillograms of <p m , 0 3 , <£ 4 , and Nl are shown in Fig. 34 for the 
case of Nl - 0.8 AT and Leg 4 unloaded. These oscillograms were taken at 
a temperature of 29°C and a repetition rate of 60 cps. Note that 
Nl E 0.8 AT is below the minimum value of Nl that enables A 4> m to be equal 
to 2<fi , as is evident from Fig. 33(a). These oscillograms were taken in 
order to explain a phenomenon associated with peaking of D in a saturable 
core, as will be discussed later. 

Our next objective is to compare measured and calculated plots 
of D versus jV I . In order to calculate these plots, we first must deter- 
mine the leg parameters. 
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FIG. 34 FLUX-DIVISION OSCILLOGRAMS OF CORE S 
WITH Nl = 0.8 AT 
Time scale = 4 /i.sec/major div.; 

0 m and 03 scale = 0.046 (V/T)/major div.; 
04 scale = 0.020 (V/T)/major div. 


C. LEG PARAMETERS 

The inelastic-switching leg parameters were measured in the 
same manner as described in Part C-4. Leg m (Fig. 22(b)] was wrapped by a 
short-circuited winding, and Legs 3 and 4 in series were tested as a 
single-leg core. The measured peak 0, 4 » is recorded versus the ap- 

plied MMF, F 34 , in Fig. 35. 

By extending the linear portion of 0 p 3 4 versus F 34 (dashed line 
in Fig. 35), it is found that F Q34 = 0.85 AT and P 3ip = 0.75 ohm /T 2 . 
Following Fig. 22(b), the leg dimensions are as follows: Z® Y = 16.845 mm, 

Z" = 5.040 mm, ZJ Y = 8.960 mm; 4 3 = 4 4 = (4^/2) = 0.665 mm 2 . Substi- 
tuting = Z® v + Z® v = 14.000 mm and 4 34 = A 3 = A 4 = 0.665 mm 2 into 

Eqs. (88) and (89), we find that // Q = 60.6 AT / m and = 15,800 ohm s/(T 2 m). 
Substituting the leg dimensions and the values of // Q and £ into Eqs. (88) 
and (89), the leg parameters in the linear region of <p p versus F are found 
to be as follows: F Qm = 1.0227 AT, F Q3 = 0.306 AT, F Q4 = 0.544 AT, 
p mp = 1.246 ohms /T 2 , p^ p = 2.083 ohms /T 2 , and p Ap = 1.172 ohms /T 2 . 

The nonlinear portion of 0 p 3 4 versus F 34 , (Sec. I-D-5-b) is ap- 
proximated by the solid line in Fig. 35. With reference to Eq. (90), p. 40, 
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FIG. 35 0 p vs. F OF LEGS 3 AND 4 OF CORE S IN SERIES 


\ J4 - 0.5 ohm/ (7’ 1 ' + 1 A v ~ 1 ) , F" 34 = 
values of l ** and 4 J4 into Eqs. 
and k = 10,520 ohms/(r v+ l A v ~ l m 2 ~ 


0.55 AT, and v = 1.3. Substituting the 
(91) and (92), we find that ff" = 39 . 3 AT/m 
v ). Hence, the leg parameters in the non- 


linear region of 4> ^ versus F are as follows: F" 0m = 0.662 AT, 

^03 = 0.198 AT, f 04 = 0.352 AT, and, in units of ohms/T 2 • 3 A 0 • 3 , k m =0.357 
^-3 “ 0.855 and A 4 = 0.405. 


These values of parameters were then employed in computing D 
versus NI for various load levels, to be discussed next. 


d. COMPUTED D VERSUS MEASURED D 


Computation of f 1 ux -d i vi s i on ratio D versus applied MMF NI was 
based on Eqs. (70) through (73). The ratio 4> r /4>, was measured ( cf . 

Part C-3) and found to be slightly larger than 0.9, and was rounded to 0. 9. 
Later, we shall examine the effect of on D. The computation of D 

was performed on the Burroughs 220 computer, using the ALGOL programing 
language . 



For a given value of N 2 L /R L and for each value of NI , the com- 
putation steps are described as follows: Initial numerical values of 

the parameters of each of the three legs are assumed to be F 0 and p p . 
These six parameter values are substituted into Eqs. (70) through (72) 
in order to solve for p t r a ns ce nde n t a 1 1 y . Switching time T[Eq. (D-14)], 
is then evaluated and substituted into E q . (D-12) in order to determine 
F 3 , which is tnen substituted into Eq . (D-4) in order to determine F m - 

The expression 



T/4> s 


( 122 ) 


is then used to determine F 4 [Eq. (122) is obtained by substituting 
= A0 4 /T * ( 4> r + 4> if )/T = (4> r ~ 4> 3f )/T t Eqs. (60) and (62), into 
Eq. (54)]. 

For each leg whose parameters v, A. and Fq are known, substitu- 
tion of F into Eqs (94) and (95) yields the first approximation for the 
leg parameters p ' and F q. These six leg parameters replace the initial 
values of leg parameters, and the same computation steps are repeated. 

This procedure is repeated several times until one of the variables, 
e.g., F 3 , has changed by less than 0.1% in the last cycle. The values 
of p and (] in the last cycle are substituted into Eq . (73) in order to 

compute Z). 

Computed plots of D versus NI For various values of are 

compared in Fig. 36 with experimental data. The accuracy of the experi- 
mental measurement is lower as NI becomes smaller, since when A0 4 becomes 
smaller, the noise level is more significant. The error in A 0 due to the 
noise level in the flux-reference (cf. Appendix F) is less than iO.l /N o 
maxwells, where N c is the number of turns of the sense winding. For the 
cases of (N*/ll ( ) = 0, 3.926, 1 0.66, and 40.0, N a on Leg 4 was equal to 2, 
20, 10, and 20, respectively. 

e , DISCUSSION 

(1) The calculated /) versus NI (Fig. 36,) is far from coin- 
ciding with the measured curves. However, in view of the 
complexity of flux division, these results are nevertheless 
e n c ou r a g i n g . The measured /) decreases as NI increases 
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(b) (Nl/R l ) = 10.66 AND 40.00 T 2 /fl 

FIG. 36 CALCULATED AND MEASURED D vs. NI WITH (N?/R. ) 
AS A PARAMETER 
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above, say, 1,3 AT, until, at a certain value of NI , D 
reaches a minimum. (No such minimum is reached in the 
no-load case for NI below 5.0 AT.) The computed D versus 
NI for the loaded case also has a minimum point that be- 
comes more pronounced as N^/R^ increases (see also Fig. 39 
later). When we first computed D versus NI , we assumed 
constant F Q and p in each leg, and the resulting curves 
were monotonies 1 ly decreasing for all values of N^/R^. 

Since the effect of the curvature of 4> p versus F was ac- 
counted for in computing the curves in Fig. 36, it may be 
concluded that the minimum of D versus NI is associated 
with the nonlinearity of cf p versus F . Since the assumption 
of constant F n and p is valid for high values of i VI, the 
effect of including the nonlinearity of <p p versus F is not 
to raise D at high values of N I , but rather to drop D at 
the low v lues of NI . We may conclude, therefore, that 
for values of NI larger than NI at the minimum point, D 
rises and approaches an asymptote as NI 00 . 

(2) The main discrepancy between the computed and measured 

curves is around the low values of NI , where the measured 
D reaches maximum. Obviously, the analysis in Sec. I-D 
does not account for the mechanism that causes D to reach 
maximum. This phenomenon is explained qualitatively as 
follows. There is a minimum value of NI (1.1 AT in this 
case), below which A <fi m < 20 r , even if T Around this 

low value of N I , flux switching follows a certain sequence. 
First, flux switching occurs essentially only in Legs m 
and 3, because //. < ff Toward the end of the flux 

switching in Leg 3, F 3 increases because Leg 3 reaches the 
wing of the static <fi(F) curve near saturation, while Leg m 
is near a demagnetized state. Such a redistribution of NI 
into higher F 3 and lower F m follows the requirement that 
$ m = $> 3 and the fact that A m > A 3 . The increase in F^ 

(=F 4 in the no-load case) increases F 4 beyond the static 
threshold value of Leg 4, and Leg 4 is then able to switch 
ine 1 as t ica 1 1 y . Leg 4, therefore, switches only after 
Leg 3 completes most of its switching. This conclusion is 
verified experimentally in Fig. 34, where NI was chosen be- 
low 1.1 AT in order to emphasize the delay in $> 4 relative 
to $> 3 . Coming back to the plot of D versus N I , as NI is 
increased beyond 1.1 AT, the pulse duration T must decrease 
in order to keep A0 = 2<f> r . Such a decrease in T will de- 
crease (because when NI is cut off, Leg 4 is still 

switching), but hardly affectA0 3 (because Leg 3 has already 
completed switching). The amount dropped in A<£. is equal to 
the amount raised in A^ , because A<£ 3 +A<£ 4 = 2(p r = constant. 
Such an increase inA</> 3 results from the increase in NI . 
This increase in A<£ 3 and decrease in A<£ 4 explain why 
D(=A0 3 /A0 4 ) increases with NI . These variations in A0 3 
and A0 4 can be seen in Figs. 33(a) and (c). In Fig. 33(d) 

the load is so high, that /f 4 does not exceed // t h , hence 
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this phenomenon does not occur. It is not understood why 
this phenomenon was not observed in Fig. 33(b). As NI in- 
creases further, it becomes high enough to switch all legs 
simultaneously, as assumed in our analysis, and D decreases 
with NI , At the transition between the region where D in- 
creases with NI and the region where D decreases with NI , 

D reaches its maximum value. Note that this mechanism is 
entirely different than the mechanism that causes D to 
reach maximum in Core N, as discussed previously. In fact, 
since all legs in Core N have the same cross-sectional area, 
the phenomenon described here cannot occur in Core N. 

(3) The measured value (0 r /0^) = 0.9 [hence, following Eq . (68), 

S = (1 - 0 . 9 ) /( 1 + 0.9) - 0.05263] was used in order to 
compute D versus N I , Fig. 36. One may wonder how sensitive 
D is to 0 r /0 s ; in other words, how crucial is the measure- 
ment of 0 r /0 s . To answer this question, D versus NI in the 
no-load case was computed for various values of 0 /0 , 
ranging from 0.86 to 0.96. The resulting plots are shown 

in Fig. 37. We conclude from these plots that, considering 
the agreement between calculated and measured D versus NI 
in general, determination of 0 r /0 # need not be within tight 
tolerances, 

(4) In computing D versus NI , Fig. 36, variation in leg param- 
eters with NI were accounted for. However, in doing so, 
the computation was much more complex than if the leg parameters 
were constant. In several cases, the accuracy of the results 
may be sacrificed somewhat in order to simplify the computa- 
tion procedure. To have an idea of how much accuracy is 
sacrificed by such an assumption, D versus NI in the no-load 
case was recalculated, assuming constant leg parameters ob- 
tained from the linear portion of B versus H . The resulting 
curve is compared in Fig. 38 with tfle measured curve and the 
curve whose calculation was based on assuming that F' Q and p' 
of each leg varies with NI . As expected, assuming constant** 
^0 and Pp in each l e g for low values of NI results in a sub- 
stantial error, but for medium and high values of N I , such 

an assumption is certainly justified. To get better agree- 
ment at the low values of NI , yet to keep the computation 
simpler, the 0^ vs. F might be piecewise linearized into 
two regions. 

(5) For future reference, the effect of the load on D versus NI 
was recalculated for a wider range and rounded values of 
/V*/tf L . These curves are shown in Fig. 39. Note how the 
minimum point of D becomes more pronounced as N*/R t increases. 

(6) The assumption of 5 1 was examined in the analysis, 

Sec. I-D-4-d; it corresponds to the assumption that in each 
l e £» P ~ Pp ~ constant. Accounting for variation in leg 
parameters with NI , versus NI for § 1 was computed from 

Eq. (80). The resulting plot is compared in Fig. 40 with com- 
puted D nl versus NI , assuming that 5 = 0.05263, i.e. , <p /<p * 0.9. 
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FIG. 37 EFFECT OF (<p,/<f> t ) ON CALCULATED D vs. Nl IN THE 
NO-LOAD CASE 
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( 7 ) 


The effect of Z 4 / i 3 on the asymptotic value of D as NI — 00 
is shown in Fig. 41 for various values of p /(fl, /N ? ) . 
Assuming that (<fi /<fi ) e 0.9, these plots were calculated 
from Eqs, (111), (112), (72), and (73). At a first glance, 
these plots look suprisingly linear, whereas in fact they 
are nonlinear. In order to understand why these plots ap- 
pear linear, let us first examine the no-load case. By 
differentiating the expression for D NL jyj^oo (113)] 

with respect to l /l^, the slope s of b fii versus 

Z 4 / Z 3 is found to be 




2 (~in S) 


1 - 5 2 


\§-y - §y/ 


(123) 


where 


*4 


(123a) 


Expressing §~ y and S y in Eq . (123) by infinite series, we 

obtain 


- In 8 


1 - 6 2 
25 


1 2 


fv (nM 3 


^ , c 


(124) 

For Z 4 / Z 3 >> 1, |y S| << 1, hence $, Eq. (124), may be 
approximated by 


~ 1 - S 2 

25(-inS) 


(125) 


Since S 2 << 1, Eq . (125) may be simplified further to the 

foil owing : 


~ 1 

28 (“ in 8) 


(126) 


It is evident from Eq. (125) or Eq . (126) why D Ni NI ^ W 
versus Z 4 / Z 3 may well be approximated by a straight line 
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FIG. 41 CALCULATED D N| . ro vs. (/ 4 // 3 ) FOR VARIOUS (a N? /R. ) VALUES 
=0.9 p 
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of a slope s which depends only on 8. In Fig. 41, 

8 - [1 - (</> r /0 # )]/[l + (0 r /0, )1 • (1 - 0.9) /(I + 0.9) - 0.05263. 

Thus, for (1 4 /1 3 ) » 1 in the no-load case, Eq. (126) 
yields s ^ 3.22, whi ch ag r ees we 1 1 with the correspond- 
ing plot in Fig. 41. For the cases where Leg 4 is 
loaded, the plot appears even more linear. The mathe- 
matical interpretation is complicated by the transcend- 
ental equation for Eq. (UD. However, it was 

shown 0 that loading Leg 4 has essentially the same 
effect on D as increasing I 4 /Z 3 of an unloaded core. 

Hence, the expression in Eq. (126) is even a better 
approximation for s than in the no-load case, and the 
plots are more linear. 

4. Summary 

Flux division in a nonsaturable core was measured (Figs. 29 and 30), 
and oscillograms of 0 waveforms (Figs. 31 and 32) were studied. These 
data are compared with previous treatments of this subject. 6,7,1 Flux 
division in a saturable core was measured (Fig. 33) and compared with 
calculated plots (Fig. 36) for four load values. Plots of D versus NI 
are also calculated, assuming switching parameters obtained from a non- 
linear <p p versus F curve (Fig. 35), under the following conditions: no 

load with various <f> r /<f> t values (Fig. 37); various load values (Fig. 39); 
and 8 -* 1 in the no-load case (Fig. 40). Also, assuming constant switching 
parameters, D versus NI in the no-load case (Fig. 38), and D versus Z 4 /l 3 
for NI — 00 (Fig. 41) were calculated. 

E. UNSETTING EFFECT 
1. Introduction 

The unsetting effect was described in Report 1 (pp. 96-99). This 
effect is a special case of a more general phenomenon, according to which 
elastic A0, A<£ e , in a leg driven far into saturation will cause an appre- 
ciable inelastic A<£ in (at least) two other legs. This phenomenon stems 
from the availability of two inelastic-switching paths, one for the rise 
and the other for the fall of A0 tf in the leg in saturation. Upon termi- 
nation of the drive pulse, the legs of the core were shown (c/. Fig. 93, 
Report 1, p. 97) to be left in a “stressed” state, i.e a residual H 
field, sustained by magnetic poles near the leg junctions, remains in 
each leg. Due to variations in air-flux linkage, a tight winding on the 
driven leg minimizes the unsetting effect. Our objective now is to 
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verify experimentally the variations in <fi and F predicted under two con- 
ditions (Fig. 93, Report 1, p. 97): a positive drive on Leg 3, and a 
negative drive on Leg 4. 

2. Experiments 

A commercial core, AMP Core No. 395813-1, was chosen for the experi- 
ment. The core dimensions are shown in Fig. 42. Three small apertures 
were ignored, and a fourth one was identified with the small aperture in 
Fig. 92, Report 1, p. 97. From the core dimensions in Fig. 42, the aver- 
age leg lengths are found to be as follows; * 14.8 mm; Z* T * 1.1 mm; 

ZJ T * 2.45 mm. The core parameters are as follows: 

H c ** 39.0 AT/m(= 0.49 oe)\ S v - 63.7 (AT/m)^sec ( “0 . 8 oe /use c); 

c 2<fi r j c maxwells; and F ^ 53 0.72 AT (around the major aperture). 

Experimental verification of the unsetting effect was performed in 
two separate experiments. The drive windings in both experiments were 
wound tightly as shown in Fig. 92 of Report 1 (p. 97) (instead of sepa- 
rate N § and N H windings, a single winding with N s turns was wound on 
Leg 3). In addition, a two-turn sense wind ing was wound on each leg. 



R A itr 


FIG. 42 DIMENSIONS OF THE CORE USED 

FOR THE UNSETTING EXPERIMENT 
AMP Incorporated, Core No. 395813-1; 
Dimensions in mm; 

Tolerance = ±0. 025mm 
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The conditions for the. two experiments were similar, except for the SET 
state. In the first experiment, the SET MMF was relatively small, so 
that the unsetting of Legs m and 3 was caused by pure elastic flux 
switching in Leg 4. In the second experiment, the SET MMF was high 
enough to cause some inelastic setting (supersetting effect) in Leg 4, 
as indicated by Point E in Fig. 93 of Report 1 (p. 97). The results of 
the first experiment are in the form of oscillograms and the results of 
the second experiment are in the form of plots. 

In the first experiment, the number of turns of the drive windings 
(Fig. 92 of Report 1, p. 97) were: N c = 7 , N g = 1, and N = 10. The 

drive currents were 2 P - sec long, and their amplitudes were: I c - 3.4, 

I = 34, and I = 124 and 204. Oscillograms of the drive currents and 
the resulting changes in <p during the CLEAR, SET, and UNSET phases are 
shown in Fig. 43, for the two values of unsetting MMF s , NI = 120 AT, 
and NI = 200 4T. These oscillograms are in agreement with Fig. 93 of 
Report 1 (p. 97). 

In the second experiment, the number of turns were N ^ = 7, N s ■ 2, 
and N = 14. The drive currents were 2 pisec long; the CLEAR and SET 
currents were kept constant (N C I C ~ 25 AT and N s I t = 6 AT), and Current I 
was varied. For each value of NI , 0 was integrated and A<£, relative to 
the CLEAR state, was measured at three different times (see Fig. 93, 

Report 1, p. 97): after I s is removed (Point E) , during the application 

of I (Point G), and after I is removed (Point J). The flux measurement 
was made by means of the f 1 ux- re f erence (c/. Appendix F). The resulting 
curves are shown in Fig. 44. These results will now be discussed. 

3. Discussion 

Let us now compare the experimental results in Fig. 44 with the 
variations in <p and F of each leg, which were predicted in Fig. 93, 

Report 1, p. 97. In general, Fig. 44 agrees with the predicted varia- 
tions of (p in each leg. There are, however a few points that should be 
discussed. In order to facilitate this discussion, let us establish the 
following nomenclature: Let P and Q stand for any point in Fig. 93, 

Report 1, p. 97. Then, &4> jP represents the flux change of Leg (j = m, 

3, 4) from the CLEAR state (Point 4, (p = ~<p rj ) to Point P; <p jP represents 
the flux level of Leg j at Point P; and A <p j{p _ Q) represents the flux change 
of Leg j from Point P to Point Q. 
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(b) N I = 200 amp TURNS 


FIG. 43 OSCILLOGRAMS OF DRIVE CURRENT AND FLUX CHANGE vs. TIME 
FOR THE UNSETTING EXPERIMENT 

Time Scale = 0.5 /xsec/major div.; Scale for I,- and l $ =1 A/major div.; 
Scale for I =5 A/m a jor div. 
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After the CLEAR MMF is removed, the core is set via Leg 3. Due to 
the supersetting effect (c/. Report 1, pp. 93-94), A0 a£ . > 2<£ r 3 ( ■= 21 
maxwells), A0 3f *= 2<£ r3 , and A0 4f > 0. The expe r iment a 1 A <p £ of the three 
legs agree with what had been anticipated except for the following: 

Since N 9 I was maintained constant (6 AT), A$£ of each leg was expected 
to be independent of JVI. Instead, A <f> mS and A</> 4£1 have a slight droop as 
NJ increases. This droop could not be explained until we found that the 
Current pulser that had been used in this experiment had a low dc leakage 
current. This current increased from 1 to 15 ma as we increased the 
number of pulser channels that were connected in parallel. Instead of 
recording new data with the dc MMF cancelled, we decided to use the data 
because the effect was small, and because the main purpose of the experi- 
ment was to measure A^j andA<£ 4 j, which, as will be shown later, are not 
influenced by the dc MMF. Furthermore, we shall interpret the data in 
considerable detail because of the interesting phenomena associated with 
the dc leakage currents in the experiments. 

The resulting dc MMF (t.e., N times the dc leakage current) of up to 
0.21 ampere turns caused Points A [Fig. 97(b), Report 1, p, 97] to be 
biased away from the CLEAR state. This dc bias pushed Point A of Leg 4 
down and Point A of Leg 3 up (in order to balance the elastic drop in 0 4 ), 
resulting in 0 3 being clamped to the static loop above -0 r3 . (The 
effect on was negligible due to the relative length of Leg m. ) The 
abscissa of Leg j (j - m, 3, 4) in Fig. 44 corresponds to the biased 
Point A, not to <fi = * SET MMF pulse, N S I S , is applied next. Because 

of the higher initial the resul t ing A0 3 ^ A _ c ^ is smaller than the 

corresponding flux change without the dc MMF. This causes and 

hence A0 b ^ j , to be smaller as well. Due to the dc MMF, relaxes to 
Point E slightly higher than without the dc MMF (causing to be lower). 

There are, therefore, two opposing effects of the dc MMF on A< ^3(4-£)- N ° 
net effect could be detected experimentally. This is evident from 
Fig. 44(b), where A</> 3 ( A _ E y is constant in NI ( NI varies with the dc 
leakage current). An additional drop in 4> mE and was caused during the 

interval between the / and I pulses by the combined dc MMF of the leakage 
current in the N{ = 14) turns on Leg 4 and the A^ c ( = 7) turns on Leg m. This 
combined dc MMF of up to 0.315 AT acted against a soft threshold around 
the major aperture. Given sufficient time (as was the case in the experi- 
ment), some of the inelastic flux that had been set into Leg 4 due to the 
supersetting effect was switched back via Leg m. 
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Compared with the case of no dc MMF, <t> XE is lower, thus 1A0 4(£ _ c) I 
is smaller because of a smaller inelastic switching. This results in 
less unsetting of Leg m, i.e., h<P m{E -j) I becomes smaller. It was found 
experimentally that A was independent of the dc MMF. This may be 
explained by the fact that the dc MMF decreases both A 4> m{A . E) and 

|a<^ ( k- j ) I ■ Durin s the fal1 time> as Current 1 dr °P s t0 the dc leaka e e 

level 0 rises from Point G to Point J, which is slightly lower than in 
Fig. 93(c) (Report 1, p. 97); hence, 0 3J is slightly higher. 

We therefore conclude that Point J of each leg is hardly affected 
by the dc leakage in the pulser. Based on this conclusion, the curves 
in Fig. 44 were plotted as follows. For each leg,A0 J (resulting from 
the CLEAR pulse) was plotted first; A0 (f _ J( was then used to plot A0 g , 
andA0 (s _ c) was used to plot A 0 C - 

The plots of &4> mJ and A 0 3j . , Figs. 44(a) and (b), show the magnitude 
of the unsetting effect. For example, Leg 3 in a SET state can be de- 
magnetized (i.e., 0 3 dropped from +0 r , to zero) by driving Leg 4 by 
NI = 128 AT (which is about 178 times around the major aperture). 

For NI = 242 AT (=336 F 0 ) , A0 jy » 5.2 maxwells, thus 75% of 2 0 r3 was unset. 

4. Fall-Time Effect 

The effect of winding position on the unset flux was discussed in 
Report 1 (p. 99). It was shown there that the tighter the drive winding 
is, the smaller is the amount of unset flux. 

The amount of unset flux is also influenced by the fall time of the 

drive current. This stems from the fact that Leg 3 is unset ine 1 as t ica 1 ly , 
and inelastic switching requires a longer time than elastic switching. If 
the fall time is too short, the (fast) elastic switching in Leg 4, being 
unable to close fully through Leg 3, would close partially through the 
surrounding air. This would result in a creation of pole distribution 
near the leg junctions. The polarity of the poles is similar to the one 
in Fig. 87(c) of Report 1 (p. 92). These poles act on Leg m in the posi- 
tive direction, and on Legs 3 and 4 in the negative direction. As a re 

suit, the excess II in Leg 3 would increase, enabling Leg 3 to switch 
faster. However, if the (soft) threshold in Leg a is exceeded, a portion 
of the elastic A0 4 is steered into Leg m in the positive direction, and 
the pole density would decrease. The net result is a decrease in the un- 
set flux in both Leg m and Leg 3. In conclusion, the unsetting effect is 
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expected to be smaller as the fall time becomes shorter (below some thresh- 
old value). This conclusion has not yet been verified experimentally. 

5. Summary 

The unsetting effect, which has been described in Report 1, is veri- 
fied experimentally (Fig. 43). After the tested core is cleared, Leg 3 
is set by a fixed SET MMF, and this is followed by applying a variable 
negative NI to Leg 4. After the removal of a relatively high NI , the 
core is left in a state commonly obtained by a relatively low NI through 
the major aperture. Additional effects due to leakage dc current in the 
driver pulser cause the set flux in Legs m and 4 to be slightly lower 
than without the leakage current. These effects are explained and taken 
into account in plotting the experimental data in Fig. 44. 

Since the unsetting of Leg 3 during the fall of the drive current I 
involves inelastic switching, the duration of the fall of I influences 
the amount of the unset flux. The shorter the fall time is, the smaller 
is the unsetting effect. 

F. CONCLUSIONS 

In discussing and evaluating the results of much of Sec. II, it has 
been assumed that the magnetic properties of a core, or a leg of a core, 
are determined by integrating the corresponding material properties over 
the x coordinate (c/. Report 1, Fig. 19, p. 18) of the leg. This assump- 
tion is invalid if the material properties at one x value in the leg are 
affected by the process of switching or the magnetic state at other x 
values in the leg. For example, if the material of a toroidal core has 
a large magnetostriction constant, then a partially switched annular 
region of the core (where |b| < B r ) , which is not free to change its 
dimensions, may produce a stress that can change the material properties 
in other regions of the core. However, we have no experimental evidence 
to support such an effect, and further investigation is needed. 

The <p{F) testing of blank ferrite disks before the ultrasonic cutting 
of cores is an easy means of detecting disks having relatively abrupt 
radial variations in material B(H) property, and of making comparisons 
between disks. Even though 0(F) testing cannot guarantee that a disk is 
entirely free of material variations, it is valuable for detecting ap- 
preciable nonuniformity. 
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Two possible modifications of the parabolic switching model were 
examined numerically in Sec. II-A-3. One accounted for P p (F) variations 
and the other for F Q {<p) variations. Analytical functions have not been 
given for p p (F) and F 0 (<£>); this remains to be done in the future. Modi* 
fication of the parabolic model is needed primarily in the si ow- s wi t ch i ng 
range. On the basis of <p(<p ) curves for a thin ring of OD/ID = 1.06, the 
result of performing these two modifications numerically was examined. 

It is concluded that some discrepancies still exist, but that very sig- 
nificant improvements do result. In some cases it may be desirable, for 
simplicity, to include only the p (F ) modification. This results in 
the same modification used in Sec. I-D-5, except for the terminology, 
i.e,, p* and F' are functions of F . 

p u 

The discussion of flux switching in a core having a re-entrant shape 
makes use of the concept of a switching line . This concept is also useful 
in discussing switching in other types of non-toroidal cores. Single- 
aperture cores with re-entrant shapes can be divided into two classes. 

In cores of the first class, the shortest closed line that can be drawn 
inside the core and around the aperture does not touch the outer edge of 
the core. In cores of the second class, this shortest closed line touches 
the outer edge of the core. The discussion on slow flux switching in the 
first class (c/. Part B-3) is not detailed; it only presents the general 
configuration of the switching lines. The discussion of the second class 
includes a proposed model, which is called the switching crossover zone. 
This model, although self consistent, is intended only as an example of 
the kind of switching that may occur. The validity of the concepts pre- 
sented in discussing re-entrant shaped cores rests upon two assumptions: 

(1) that zone boundaries actually exist and have approximately 
the properties described in Report ], pp . 119-128, and 

(2) that the static B(H) loop of the ferrite material has 
vertical sides. It is hoped that the basic concept of 
the crossover zone model will be useful in studying re- 
entrant shaped cores made of ferrite having more realistic 
static B(H) properties. 

The section on design considerations was specifically concerned with 
the design of two cores to be used for experimental studies of flux 
switching in multipath cores. The basic goal was to design two cores 
that would exhibit as few superfluous phenomena as possible. These 
considerations may not be applicable to the design of cores for optimum 
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performance in a practical circuit. The final designs are not necessarily 
the best that could be achieved. For example, in determining l I q and 
Z av of each leg of Core N after the core was cut, the problem of the re- 
entrant core shape was first realized. Thus, the design of Core N could 
be improved by making Leg 2 straight instead of slightly curved. 

The hyperbolic B(H) model, which was used to obtain the static 0(F) 
equations in terms of core dimensions and material parameters, is experi- 
mentally checked for a core made of Telemeter Magnetic T-5 material [ cf . 
Fig. 28(c)], The disagreement in 0 for “5// . Z < F < + 5H . Z is reason- 

t n o t h o 

ably small, when compared to <p r . However, the disagreement in 0 from 
~0 r up to 0 corresponding to F c // th Z fl is about 50% of the measured A0 
[cf. Fig. 28(d)], This means that the hyperbol ic B ^ (//) function, which 
fits the purely elastic switching region < H < 0, does not describe 
B(H) properly in the region 0 < H < // 1 h . It is evident that this region 
is characterized not only by elastic switching (as assumed by the model), 
but also by inelastic switching (of the same order of magnitude as the 
elastic switching). If it becomes important to describe B(H) accurately 
in the region 0 < H < then an additional function will have to be 

introduced. Evaluation of the hyperbolic B(H) model in describing the 
static B(H) properties of other materials has not been made. 

The analysis of flux division was based on two factors: variation 

of p with 0, and the nonlinearity of p and F Q at low values of F in each 
leg. From the comparison between calculated and measured D versus NI , 
it seems that these two factors are major contributors to the phenomenon 
of flux division. However, there are other factors (not accounted for 
in our analysis) that introduce errors to the calculated plots, especially 
when slow switching is involved in Leg 4 or all legs. These factors are: 

(1) A small amount of inelastic flux change at a field lower 
than the threshold value, 

(2) The limitations of the parabolic model of p in describing0 
at low F, even if the nonlinearity of B versus H is ac- 
counted for, and 

(3) Elastic flux change during the rise time of NI may differ 
from the elastic flux change during the fall time of NI . 

Inclusion of these factors in the analysis is bound to yield better re- 
sults, but this involves establishment of additional switching models. 
Rather than improving the analysis of flux division in a saturable core, 
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we feel that our next step should be to establish switching models for 
a soft state, so that we can analyze flux division in a nonsaturable core. 

Experimental verification was given which clearly substantiated the 
existence of the unsetting effect and gave the magnitude of the effect 
for a particular core. Calculations of this magnitude, based on geometry 
and material properties, remains to be done. 
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APPENDIX A 

CALCULATION OF STATIC 4>(F) CURVE 
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APPENDIX >1 


CALCULATION OF STATIC 0(F) CURVE 


The objective is to c a 1 cul at e a s t at ic 0(F) curve of a cons t ant ** width 
leg; the calculation is based on the suggested hyperbolic models for B(H) 
[Eqs. (10) and (11)] in saturation and between saturation levels. This 
calculation is done by substituting Eqs. (2) and (3) into Eqs. (10) and 
(11) in order to derive B{F t x), which will then be integrated from x » 0 
to x = w [Eq . ( 7 ) ] . 


For H < tf th » # is equal to 


Vo F 


_D — ( D _ D \ - 

r s 1 r F ~ H (l . + ax) l . + ax 


(A- 1 ) 


For H ff th . B is equal to 


f ~ H (1, + ax) p. 0 F 

B » ~B + (B + B ) + 

“ r ss r F - H ( 1 . + ax) l. + a x 


( A- 2 ) 


In deriving the expressions for 0(F), three regions of F are dis- 
tinguished: F < l ( H th <F< l o H th> and I 0 ff t h < F. 


1- T i 


For this region of F, H < // th throughout the leg cross section, and 
only Eq. ( A- 1 ) is valid. By substituting Eq. (A-l) into Eq. (7) and 
integrating from x = 0 to x = w , we obtain the following expression: 


0 = 


Fh 


'B ~ B ( F ~ H l \ fl\ 

-*T "{thcO ' Mi;, 


- B _ hw 


(A- 3) 
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2 . 


III. < F < l II . 

i t k — — o t li 


The leg cross section is divided into two sections: In the inner 

section, H > // . , and in the outer section, II < H Let x be the 

distance from the inner edge to the boundary between the two sections. 
At x = x ft H = // th , hence Eqs. (2) and (3) yield 


— (~~ ~ i ) 
a Uth *' 


( A- 4 ) 


In the section where 0 x _< , Eq. (A-2) is valid; in the section 

where x f x w , Eq. (A-l) is valid. Therefore, Eq. (7) becomes 


4> = h[J B (F,x )dx + J BAF,x)dx ] 

0 »* 


( A- 5 ) 


After substituting Eqs. (A-l) and (A-2) into Eq. (A-5) and integrating, 
we get the following expression: 


* ■ 

a 


B - B 

st r 

H 


In 


F - H l 

a 0 


I H «' 

F I- - 
B *ki 


+ In 

H 2 


H 


F 1-; 


H 


t h/ 


F - H l. 

n i 


(B + B)H 

s s r a 

+ 1 + 


( ^ O \ 


"(B„ +B r W 


► - ft 

+ B w 

\ l i) 


L 7/ n a ' J 


(A-6) 


3- 1 F 

Here H > H t h throughout the leg cross section, and only Eq. (A-2) 
is valid. After Eq. (A-2) is substituted into Eq. (7), and integration 
is performed from x * 0 to x - u> , the following expression is obtained: 


0 


fh 

a 


H 2 


In 


F-HJ,, 




/n 


~(B ss +B r )H 

j — ] 


- B 

\U 




( A- 7 ) 
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APPENDIX B 


RELUCTANCE OF A VARIABLE -WIDTH LEG OF LINEAR MATERIAL 


Let a point in space be described by rectangular coordinates ( x,y t z ), 
and by orthogonal curvilinear coordinates (u^u^u^). There is a func- 
tional relation between x, y , and z and tij, u 2> and u j. It can be shown 5 
that an element of distance is expressed by 


where 



(B-l) 


(B-2) 


Since we have assumed a leg to have a constant thickness h , our problem 


is two dimensional 

, and 

Eqs. (B-l) 

and (B-2) are reduced to 



ds ~ 

J BlW”! * 

S 2 2 du l 

(B-3) 

where 

6ll = 

( 


(B- 4 ) 


\ Buj / 

VBu,/ 

and 

S 2 2 


£)’ • 

(B-5) 

Now consider 

a var i 

ab 1 e - wi d t h 

leg whose boundaries follow 

the direc 


tions of some orthogonal curvilinear coordinates Uj and u 2 - In Fig. B-l, 

the edges of such a leg are defined by u 2 * u 2a and u 2 = u 2fc , an< ^ 

ends of the leg are defined by Uj = u la and u } = u 1 , where u 2a , u 2b * 

u . , and u , , are all constants, 

la’ 16 
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EQUI-F LINE 



FIG. B-l CURVILINEAR LEG 


The leg is now divided into curvilinear squares of length vg ^ i ^ u j 
and width ' / g 22 du 2 that, together with thickness h , define leg elements. 
The length and width of each curvilinear square vary in such a way that 
the MMF drop along the leg element, HV g { } du j, and the flux through it, 
Bh^ g 2 2 du 2 > are identical for all the leg elements. By plotting these 
curvilinear squares, one gets the equi-F lines along the u 2 direction 
and the flow lines along the u } direction. The reluctance of every leg 
element defined by a curvilinear square is the same, and following 
Eqs, (26) and (27), i s equa 1 to 


Hv'g U du. 

AR = _ 

Bhvg 22 du 2 

The reluctance of i>R of a larger leg 
equi-F lines, is the inverse of the 


1 

‘ = — • (H-6) 

M 1 Z 2 2 du 2 

element, defined by two adjacent 
integral of the inverse ofAR, i.e., 


SR 


du 


) 



(B-l) 


The reluctance of the whole leg is obtained by integrating SR along the 
leg; hence, 
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/ n 


R = — 

fih 


da , 


2 fe 


p 2 2 


dli 


u , u 


« 1 1 


0 2d 


Let us now apply Eq. (B-8) to three examples. 

Example 1 — A toroid of inner radius r. and outer radius r o . 

Here, u t - 0 {a ^ a - 0; u j ^ = 277); Uj ” r ( u 2a ~ r t * u 2 6 * r o^ 
x = r cos 6; and y = r sin 0. Following Eqs . (B-4) and (B-5), gjj 
and g 2 2 = 1 . Thus , 


R - 


1 

jjih 



R = — 


l 

jjih 


2v 


l n 


7 ) 


Example 2 — A tapered leg (see Fig. 3) of length r fc - r a and angle / 


Here, u. = r (u. - r ; u 

* l I a a 


1 b 


) ; u n ~ 0 ( u « = 0; n 


- r cos 6\ y = r sin 9\ g^ = ]; g 


2 2 


k 2 ~ v 2 a 

r 2 . Thus . 


2 fc 


/*>; 


l 

fih 



(B-8) 


( B- 9 ) 
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E 


Example 3 — A t rans f luxor - 1 i ke leg, Fig. B- 2 , defined by angle a and 
eccentric circles of radii r a and and eccentricity e. It can be 
shown 5 that the circles can be described by bicylindrical coordinates 
(77,6?) that are related to the rectangular coordinates according to 

a sinh 77 

x = , (B-l 

cosh r} - cos 6 


a sin 9 


cosh 77 - cos 6 


where a is a scale factor which is related to r , r. and e as follows 

a 6 


2 e 2 ( r ? + r?) + (i 


FIG. B-2 TRANSFLUXOR LEG 




A flow line follows a circle of radius 


a 

|sinh 17 1 

with a center on the x axis at a distance 

d ~ a coth 77 


(B- 14) 


( B- 15) 


from the origin. By varying 77 from tj * to (r) a > T)^) , the circles ex- 
pand and move to the right, i.e., the radii increase from r a to and 
the centers move from x - d a to x = 

A constant-F line, perpendicular to the flow lines, follows a circle 
of radius 


| s in 8 | 

with center on the y axis at a distance 

d' = a cot 8 


(B- 16) 


(B- 17) 


from the origin. The leg ends are thus sections of circles defined by 
letting 8 = ±0£ in Eqs. ( B- 16) and (B-17). 

In order to calculate R, let us divide the leg into two symmetrical 
half legs, one above the x axis and one below it. The width of each half 
leg varies as 8 varies from ±77 to iOt . The total R is twice that of one 
of the half legs. For the top half leg, u l - 8 (« l3 - ; a 1 = Of ) ; 

u 2 s 77 (u 2q = T ?a ; u 2fc * Tj h ) . Following Eqs. (B-4), (B-5), (B-ll), and 

(B- 12) , 


^ 22 (cosh 77 - cos d) 


The total R is then 


R = 2 


fih 


dd 


J 4 dr, 

77 TJ 


(B- 18) 
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or 


1 2(7 t - a) 
M V a ~ V t 


(B- 19) 


Substituting Eq.‘(B-14) into Eq. (B-19), 


R = 



(B-20) 


where a is expressed in terms of r fl , and eccentricity e in Eq* (B-13). 

As a final note to the problem of finding R in general, suppose 
that the leg boundaries are not amenable to analytical expression. In 
a case like this, we shall have to resort to numerical or graphical ap- 
proximations. 4,5 The smaller the curvilinear squares are, the more 
accurate is the approximation, but this makes the process more tedious. 

It is better to start with large curvilinear squares so as to have a 
feel for the general trend, and then reduce their size step by step. 

Suppose we end up with rij x n 2 curvilinear squares, each with the 
same length ratioA lj/A l 2 (subscript 1 corresponds to the direction of 
flow lines along the leg; subscript 2 corresponds to the direction of 
equi-F lines across the leg). The total R is, then, 


R 


1 A i i n i 
/J-h A l 2 n 2 


(B-21) 
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APPENDIX C 


The expression 


The expression for p 


P 


where c = r./r and 

o a 


P 


The integrand in the 



CALCULATION OF p t 


for P(uj) in general [Eq. (47)] is 



t was shown to be [cf. Eq . (33)] 



u = 4>/4> f . Substituting Eq. (C-2) into Eq. 



denominator of Eq. (C-3) can be written as 


= In (c + u) + In (c _ it ) ~ In (1 + u) 


- In (1 - u) 


(C-l) 

(C-2) 

(C-l), 

(C-3) 

(C-4) 
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Hence , 



da 


( c + u) In (c + u) “ (c ” u) In (c ~ u) 

” (1 + u) In (1 + u) + (l~u) In (1 ~ u) ( 

(C-5) 


or, 




' (c + u) C + M (l - u) 1 ~ u ~ 
Ac - u) <_,, (l + «i) I + u - 


(C- 6 ) 


After substituting the limits of integration u 0 and into Eq. (C-6), 
Eq. (C-3) gives 


P t ■ 2C p ^ 


u f~ u o 


In 


„ C +U 


1 -U 


(c + ) ^(1 ~ U^) ^(c ” Uq) °(1 + u Q ) 


l + u, 


■ (C-7) 


1+CI ■ 


c - u^) f (l + U f ) He + u 0 ) °(1 - a Q ) 


1“U / 


Substitution of u 0 = ~u r into Eq. (C-7) results in Eq. (50). 
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APPENDIX D 


CALCULATION OF FLUX -DIVISION PARAMETERS 


The six basic equations for calculating flux-division ratio D are 
Eqs . (52) through (57). Two of these equations have already been employed 

in the text. Equation (52) was used in order to show that c 3 
[ cf . Eq. (62)]. Equation (56) was used to derive Eqs. (65) and (66) 
which, when combined, become 



(D-l) 


f 1 

Equation (D-l) includes two unknowns, T and J F,dt; hence two more equa- 

0 J 

tions are needed. These two equations will be derived from four of the 
basic equations that have not been used yet: Eqs. (54) and (57) will 

lead to one equation, and Eqs. (53) and (55) will lead to the other. 

For the sake of convenience, let us rewrite Eqs. (54), (57), (53), 
and (55) in this order. 


and 



(D-2) 


(D-3) 


NI = F + F, 

■ 3 


(D-4) 
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Substituting the expression for F 4 [Eq. (D-2)] into Eq. (D-3) results 


i n 


*4 


F - F 

r 3 r 04 


Nf 


P Kp /<£, ) 2 ] 


Hence 


(D-6) 


4/ 

{ 


L P ip l 1 - (0 4 /0,) 2 ] 


1 

h J 


d(p * S_ (F 3 - F 04 )dt 


(D- 7) 


After integration and substitution of . = -0 [Eq. (62)] , Eq. (D-7) 
yields ' 


4> 

P 


~ [ tanh " 1 ('^)" tanh " 1 ('^ / ) *T l (0 '"V " ■ (D-8) 


, n l 


Equation (D-8) is the second of the equations needed to solve for 4> /<p 

The third equation will be derived as follows. 

Following Eq. (D-4), 



NIT ~ J T F dt 

0 " 


(D-9) 


Following Eq. (D-5), 



(D- 10) 
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hence 


l F . d < ■ F o. 


2*. 

T + tanh 1 l — 1 

■ n ' d> / 


(D- 11) 


Substituting Eq. (D-ll) into Eq. (D-9), 


T 2 
l m t..h ‘(-J 


m p 




(D-12) 


Equations (D-l), (D-8), and (D-12) will now be employed to solve 

for 7 [which is equal to 2 tanh -1 (<P 3f /<P s ) ]■ Substituting the expression 
for tanh -1 (<P 3f /4>,), Eq. (D-l), into Eq. (D-8), We get 


V>. 


'*r\ N l 


3 p 


— 1 tanh -1 1 — I + — (<p ~ 4> 3f ) = U + 

p» W R i \ 




(D-13) 

Substitution of Eq. (D-12) into Eq. (D-13) yields an expression for T : 


2 <£> tanh 


-1 



1 1 


3 /> 


"l 


— + — 1 + — +— W>, 

Pup \ P» } fl * 


L P *p 






— . (D-14) 


m - f .- ) (‘ t ^f)' 

r ip 


^04 ^03 


ip 


Subtracting from both sides of Eq. (D-12), and substituting the 


resulting expression for (J F-dt ~ F Q3 T) into Eq. (D-l) results in 



(D- 15) 


Substitution of the expression for T [Eq. (D-14)] into Eq. (D-15) results 
in the expression 
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So * - 


t anh 


- 1 



a 


3 


(D-16) 


in which a p a 2> and a 3 depend on core parameters and NI as follows: 


N *p ip 

a i ’ P 3p m ~ F o. ~ ^03> ■ (D- 17) 

n L 

'•P / 

(D- 18) 

and 

a 3 = {NI ~ F o.^P 3p + P ip ) ~ F 0i P ip ~ F 03 p 3p . (D-19) 


Equation (D-16) is transcendental, and <f > 3 ^/ < ^ > , is solved for numerically 
after the coefficients fl j , Oj* and a 3 have been evaluated. Following 
Eq. (68), 


t anh” 1 I 



<t>r 

1 + 

1 / <t>. 

~ in, 

2 \ <P r 

1 


In S 


(D-20) 


Substituting Eq. (D-20) into Eq. (D-16) and multiplying both sides by 2, 


4> r 4> 


2a. — 

V <*> 


3 / 


0. 


a 2 In S 


7 


(D-21) 


or, 


7 ** p ~ q In B 


(D-22) 
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in which 


and 



(D-23) 



(D- 24) 


After substituting the expressions for flj, a 2 , and o 3 [Eqs. (D-17) through 
( D- 19)], Eqs. (D-22), (D-23), and (D-24) are identical with Eqs. (69), 
(70), and (71), respectively. 
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APPENDIX E 


SECOND -ORDER INTEGRATOR 


The measurement of flux changes or observation of <p(t) waveforms 
requires the use of an integrator. The commonly used R-C integrator has 
the problem of integrator droop. One way of overcoming this problem is 
to add a signal that is approximately equal and opposite to the error of 
the R-C integrator. 2,12 This error-correcting signal can be obtained by 
integrating the output of the R-C integrator with another R-C integrator, 



FIG. E-l SECOND ORDER R-C INTEGRATION 


as shown in Fig. E-l. The wide-band preamplifier has a high input im- 
pedance (* 1 megohm) and a low output impedance (less than 100 ohms). 
The outputs Ke i and e 2 can be added algebraically in a dual-trace 
oscilloscope preamplifier. Assuming that Capacitors Cj and C 2 are 
initially discharged, the following differential equations illustrate 
the principle of operation: 


Ke 


K 

fljC, 


edt ~ 


K 

filC, 


e j dt 


and 


K 


^ 2 *^ 2 


€ ,dt ” 


^ 2^2 


e 2 dt 


(E-l) 


(E-2) 
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By making fljCj = R 2 C 2 and adding Eqs. (E-l) and (E-2), we obtain 


Ke 


K 


1 

ecft _ 

«iCi 


e 2 d t 


( E- 3 ) 


Note that the amplified error of the first integration, ~(K/R jC j)J e jdt , 
has been cancelled by the output of second integration, except for a 
second-order error, - J(l//? 1 C 1 )e 2 dt. The main signal, which contains the 
first-order error, is Ke j , and the error- correct ing signal is e 2 . A com- 
parison of experimental integrals for a first-order and a second-order 
integrator is shown in Fig. E-2. The value of RC for this figure is 

300 /isec. 

Further correction can be obtained by adding more stages, providing 
a means is available for adding three or more signals that have a common 
ground . 

Greater accuracy may also be achieved by slightly overcorrecting the 
error with a second-order integrator. This can be accomplished either 
by making R 2 C 2 < R X C x or by having unequal gains in the two channels of 
the preamplifier in which addition is accomplished. The second-order 
integral in Fig. E-2 was not overcorrected. 



INPUT, e 


OUTPUT 


— SECOND -ORDER, Ke , + e 2 

- FIRST-ORDER , Ke, 


TIME — 

RP- 3696- 201 


FIG. E-2 OSCILLOGRAM COMPARING FIRST-ORDER AND SECOND-ORDER INTEGRATION 
RC = 300 fisec; Time Scale = 50 /isec/major div. 
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APPENDIX F 


FLUX REFERENCE* 


The flux reference is an instrument for measuring either remanent or 
total magnetic flux switched by a current pulse. The use of the flux 


reference 

overcomes the problems of 


(1) 

Not having a fixed flux standard 

(2) 

Integrator droop 


(3) 

Amplifier gain variations and 

nonl inear i t ies 

(4) 

Oscilloscope cat hode- ray- tube 

di s t ort i ons 

(5) 

The tedious job of accurately 
on the oscilloscope. 

reading signal deflections 


The flux reference consists basically of a standard core which is 
switched (from negative to positive remanence) simultaneously with the 
core to be measured, and a special integrator that integrates and sums 
(algebraically) JV0 of the standard core and the unknown N<p. The standard 
core has many separate sense windings that can be connected in series in 
order to cancel out the EMF due to the flux switching in the tested core. 
The circuit is shown in Fig. F-l(a). The component values given in this 
figure are the values for the flux reference that was built in this 
laboratory, and are included only as an example. Potentiometer R^ gives 
a continuous fine variation of the flux reference output. 

The R-C integrator consists of two resistors, R^ and R 2 , feeding 
into the same capacitor C. If Capacitor C is discharged at t * 0, then 
at t * t 0 the output e c is given by 


e 


c 



8 t d 



e c dt 


( F- 1 ) 


* 

Developed by E. K. Van De Riet of Stanford Research Institute. 
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INPUT FROM SENSE 



(a) CIRCUIT 



TYPICAL 

SENSE 

TURN 


(b) WIRING OF FLUX- REFERENCE CORES 


FIG. F-l FLUX REFERENCE 
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where e 


s t d 


is the output EMF from the standard core, e is the unknown 


EMF whose time integral is to be measured, and P is the total resistance 
on each side of C. The quantity e t d is given by 


8 t d 


W + P 4>^td 


(F-2) 


where N is the number of turns connected on the standard core, P 4 is the 
relative setting (0 < < 1) of Potentiometer P 4> and <& atA is of the 

standard core. The value of P is 


P ’ P i + fl 4 P 4 (l - P 4 ) = r 2 + fi 3 + fl s P 5 (i - P 5 ) 


( F- 3 ) 


. . 'o f o 

If a positive J e dt is to be measured, then J e A dt must be negative, 

0 u 0 8 1 d 
and /V and P 4 are adjusted so that e c = 0 at t 3 t Q . Under this condition, 

Eq. (F-l) reduces to 


, 0 r* 0 r 0 

J e eft * “J e ,dt + 2 J e r dt 
i) o 8td o c 


(F-4) 


r u 

The quantity j is obtained from Eq. (F-2), i.e., 


€ * t 


(/V + P 4 )A0 


s t d 


(F-5) 


whereA0 atJ is known from having calibrated the standard core. The term 
, t o 

2J e c dt is the error in the measurement. Note that the error is pro- 
portional to the net area under the e^it) curve between the times t - 0 
and t s t Q [c/. the flux sum waveform in Fig. F-2]. This net area can 
be reduced to zero by adjusting the timing of the f lux- re f erence SET 
pulse relative to the SET pulse of the core to be measured, thereby 
eliminating the integrator droop problem. The area under the e c (t) 
curve for t > is of no concern. The effective time constant of this 
integrator is PC/2, as shown by the coefficient of the last term in 
Eq. (F-l). Potentiometer /? 5 is ganged to Potentiometer /? 4 (P 5 = /? 4 ; 

P 5 = P 4 ) in order to equalize the time constant variations (as P 4 varies) 
for both inputs to the P-C integrator. The maximum t ime - c ons t an t varia- 
tion [c/. Eq. ( F - 3 ) ] withP 4 is about 10.2%. However, as long as both 
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FIG. F-2 FLUX AND CURRENT WAVEFORMS IN MEASURING INELASTIC &d> 
BY FLUX REFERENCE 


time constants vary together, no error is introduced. The time constant 
of the two inputs are made exactly equal by switching and S 2 to their 
lower position and adjusting /? 3 to give zero output [cf. Eq. ( F- 3 )] 
Second-order integration could be used, thereby reducing the accuracy 
to which the plus and minus areas of the e c (t) curve must be made equal. 
Even without second-order integration, the areas need only be crudely 
nulled. This can easily be judged by looking at 4>(t) on the oscilloscope. 

The standard core, which is a 30-50 mil ferrite memory core, is 
connected to another identical core (A0 £ - cance 1 1 a t ion core) so as to 
cancel the elastic flux changes of the standard core during SET time 
[see Fig. F-l(a)]. The edges of these cores were rounded so that wiring 
short-circuits could be avoided. After the 42 turns were wound, the 
cores were potted for permanency of wire placement, and placed in a 
small temperature-controlled oven to remove problems due to temperature 
variations. The temperature was set well above ambient, so that auto- 
matic temperature control by heating only could be employed. The oven 
heater operates from a 28 v dc source. 
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The R-C network in the CLEAR winding serves to increase the CLEAR 
pulse fall time to about 1 /zsec in order to eliminate the possibility 
of unsetting due to ringing. 

The CLEAR and SET pulses required for the flux reference can be 
generated either externally or internally, so long as the pulse timing 
can be properly adjusted. An internal pulse source, of course, makes 
for a simpler experimental setup and reduces the possibility of errors 
for routine experimental measurements. 

The units of the flux reference are not standard units of webers 
or maxwells, but need to be calibrated in some such standard units. 

This calibration process must be very carefully done, unless only rela- 
tive flux values are to be measured. The instrument used in this labo- 
ratory was calibrated by carefully photographing the <^>(f) waveform of 
the flux reference and measuring its area. The resulting calibration 
factor was: one unit = 3.64 maxwells (or 0.0364 v /zs ec per turn) ±1%. 

The noise level typically limits the precision to within ±0.1 maxwell. 
The maximum capacity is 110 maxwells. 

The actual use of the flux reference to measure the flux switched 
by a particular drive pulse consists of three steps. First, the flux 
reference is adjusted to obtain a null on the oscilloscope. Second, the 
time delay of the flux reference drive pulse is adjusted until the flux 
sum waveform has approximately equal positive and negative areas. (In 
some cases, the first step may need to be repeated. ) Third, the flux 
value is read from the dials of the flux reference. 
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APPENDIX G 


TECHNIQUE FOR OBTAINING 4>(<t>) 
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APPENDIX G 


TECHNIQUE FOR OBTAINING 4>(<p) 


Automatic tracing of <p(<p) can be accomplished by feeding the <p 
signal to the vertical axis of , an oscilloscope and the integral of (pit) 
on the horizontal axis. The integration of <£(t) to obtain (pit) for the 
horizontal signal can satisfactorily be accomplished by the use of 
second-order integration (see Appendix E). The circuit is shown in 
Fig. G-l(a). The <p preamplifier is used to introduce a time delay that 
is equal to the time delay obtained in the <p preamplifier. Figure G-l(b) 
shows the timing of the pulses, and Fig. G- 1(c) shows the resulting <£>(0); 
the numbers correspond to the numbers on the (pit ) curve in Fig. G-l(b). 
The line segments 4-5, 7-8, 11-12, and 14-1 in Fig. G-l(c) are a result 
of integrator droop. They are exaggerated in the figure for clarity. 
Since the CRT beam spends a relatively large amount of time on the seg- 
ment 14-1, this portion of the trace is overly bright. A Z-axis pulse 
is used to eliminate this and other undesirable portions of the trace. 

By the use of multiple exposures, the function F{<p) can be superimposed 
upon an oscillogram of <p(<p) , thereby indicating the part of <pi<p) for 
which F{<p) is constant. 
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APPENDIX U 


TEST PROCEDURES FOR DETERMINATION OF LEG PARAMETERS 


Step-by-step laboratory test procedures for determination of leg 
parameters, as discussed in Sec. II-C, are described as follows. 

1. STATIC <p(F) EXPERIMENT 

The circuit diagram is shown in Fig. H-l. Inductor L is an air core 
inductor, which serves to block any current that results from N s <p during 
the CLEAR time. The flux change can be measured either by means of the 


CLOSED MAGNETIC PATH 



CLEAR PULSE 

FIG. H-l STATIC 0(F) EXPERIMENT 


flux reference (c/. Appendix F) or by means of a second-order integrator 
(c/. Appendix E). The test procedure is outlined as follows. 

a. Mixed Inelastic and Elastic Flux Switching 

(1) Select a closed magnetic path in the core to be tested 
(c/. Sec. II-C-2), preferably composed of uniform cross- 
sectional area. 

(2) Block all other paths by means of a short-circuited 
winding ( cf . Sec. II-C-2). 
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(3) Connect the circuit as shown in Fig. H- 1 , and adjust 
the repetition rate of the CLEAR pulse to about 60 cps. 

Be sure that the dc SET winding has a low enough re- 
sistance so that the SET winding would not heat the 
core for large values of I, , 

d c 

(4) Establish an I ^ that is sufficient to give a repro- 
ducible remanent state by making the pulse duration long 
enough (e.g., 10 ^sec) and the amplitude high enough. If 
the A 0 observed for a value of N s I dc slightly beyond 
threshold changes as N C I^ is increased (say, by a factor 
of 2), then the initial value of N C I C is insufficient. 

(5) Record /VA0versus tf ? I dc# adjusting N C I C at each point 

such that the quantity C I c ~ ^1^) remains constant 

(c/. bigs. 17 and 18). MeasureA</> at its peak value durinc 
the CLEAR pulse. 

(6) Be sure to take a point for /V,7 dc = (1/2 )N C I C so that 
the location of 0 = 0 on the graph can be determined, as 
will be shown later. 


6. Elastic Flux Switching 

(1) Remove I , . 

d c 

(2) Make equal to the value used in Part a above to 

clear the core initially. 

(3) Record /VA 0 versus N C I C for N C I C varying between zero and 

the value used in Part a above. Measure /VA0at its peak 
value during the CLEAR pulse. Record )VA0and N I as 
negative quantities. C C 


c. Data Processing 


(1) Shift the origin by subtracting from each AA0value 
(assumed positive) of Part a one half of the /VA0value 
obtained for NJ dc = (1/2 )N C I C . The resulting differ- 
e n c e s arc the corresponding N<p values . 

(2) Divide these values of A0 by N and multiply by the cali- 
bration factor of the flux reference or the second-order 
integrator to obtain <fi in conventional units (e.g. 
maxwells or v /isec per turn). 

(3) Plot these 0 values versus F, where F is N,I d of Part a(5). 
This curve will fall in the first and fourth quadrants of 
the graph. 

(4) Determine 4> r from this graph by reading |<£| for F = 0. 
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(5) Divide each of the AA 0 values of Part 6(3) by N and 
multiply by the calibration factor of the flux refer- 
ence, or the second-order integrator. 

(6) Subtract from each of these negative A<pvalues the 
value of <p r determined in Part c(4). 

(7) Plot these negative <p values versus F, where F is N C I C 

of Part 6(3). This curve will fall in the third quadrant 
of the graph. 

Determination of the static <p{F) curve is now completed. If neces- 
sary, a static <p{F) loop may be drawn by adding a curve that is symmetrical 
with respect to the origin. 

2. 0 p (F) EXPERIMENT 

The circuit diagram is shown in Fig. H-2. The mercury choppers and 
the calibrated dc voltage references in this figure serve to measure the 
SET pulse amplitude and the peak switching voltage Ncp conveniently and 
accurately. The choppers superimpose on the oscilloscope both the signal 
(SET pulse or N<fi) and the dc reference voltage. The reference voltage 
can then be adjusted to be equal to the pulse amplitude or the peak N<p, 

CLOSED MAGNETIC PATH 



FIG. H-2 <£ p (F) EXPERIMENT 
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and its value read from the dials of the voltage reference. This pro- 
cedure eliminates such problems as cathode - ray- 1 ube distortion, time 
variations of the base line, parallax with the graticule, time varia- 
tions of the oscilloscope vertical amplifier gains, and inaccuracy in 
reading the values from the graticule. Care must be taken to avoid 
distortion of the signal in passing through the chopper. The calibra- 
tion resistor /? must have a very low L/B ratio and an accurately known 
resistance value. The chopper technique described here can also be used 
to measure the CLEAR pulse amplitude in this experiment and both the 
SET and CLEAR pulse amplitudes in the static 0(F) experiment. The test 
procedure is as follows. 

(1) Select a closed magnetic path in the core ( cf . Sec. II-C-2). 

(2) Block all other paths (cf. Sec. II-C-2). 

(3) Connect the circuit as shown in Fig. H-2. 

(4) The rectangular SET pulse should have a short rise time, 
so that most of the switching occurs while I s is constant. 

The SET pulse duration is not critical, so long as it is 
sufficient to allow 0 to reach its peak value. 

(5) Establish an N £ I c amplitude that is sufficient to give a 
reproducible remanent state. If the0(t) observed for a 
given value of N f T j is changed significantly as ACI is 
doubled, then the initial value of N C I c is insufficient. 

(6) Record N$> p versus N t I s . Do not take 0 readings for /V I 
larger than about 80% of the value for which 0 begins'to 
occur during the rising portion of the SET pulse. 

(7) Divide the /V0 p values by N and plot 0 versus F, where F 

is N I . p 
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Air flux, 6 

Analysis of flux division, see 
Flux division analysis 

Average values of p, 17-24 
p t 17-19 
17 ~ 18 

p(<t> f ), 19-22 
p(<t> f ) plots, 21-22 

22-23, 121- 122 
plots, 23 


BUT) curve, static, see Static B(H) 
curve, model of 

B(H) loop, static, 1 

Blocking flux switching in leg, 71-72 


Coercive field, 5 

Computation procedure of flux division, 
45, 87 

Const ant- wi dth leg: 
flux in, 2 
geometry, 2 
p, of, 8, 14-16, 39-40 

p(<t> f ), of, 19-22 

Contour of leg, 60 
Core: 

machining, 60 
nonsaturable, 24-25, 59 
of re-entrant shape, 61, 64 
saturable, 24-25, 59 

Core, tested: 

MAD, static <t>(F) , 53-54 

MAD, unsetting effect, 96-101 

nonsaturable, flux division, 61-63, 

69, 76-82 

nonsaturable, static <p(F ) , 72-75 
saturable, flux division, 61-62, 69, 
83-90 

saturable, Legs 3 and 4 of, d> (F) , 76, 
85-86 , V 

thin ring, 4>{<t >) , 56-57 
thin ring, 4> p {F ) , 56, 58 
Curvilinear squares, 12 


Design considerations of experimental 
cores, 59-62 
junction problem, 60 
leg contour, 60 
leg cross section, 59-60 


Dimensions, legs: 

effect on flux division, 45-47 

effect on switching parameters, 37, 39-41 

of nonsaturable core, 59 

of saturable core, 59 

Disks, testing of, 53-55 

Division of flux, see Flux division 

Drive, effect on leg dimensions, 61, 63 


Electrical circuit analogue for flux 
division, 27 

Experiment of: 

flux division, nonsaturable core, 76-80 
flux division, saturable core, 83-88 
static 4>(F), 72-76, 149- 151. 

unsetting effect, 95-101 
<£ p (F), 73, 76, 151- 152 


Fall time, effect on: 
flux division, 101-102 
flux- reference measurement, 141 

Flux: 

effect on p, 14-16 

final , effect on p and p ^ , 20-23 

in constant- width leg, 2 

vs, MMF, static, 2-3, 5-6, 109-110 

Flux density: 

nonuniformity of, across leg, 8-9 
vs. H, inel astic, time rate of. 8, 37-39 
vs, H, static, see Static B(H) curve, 
model of 

Flux division: 
in Laddie, 24 
in loaded core, 25-26 
in MAD, 25 

in nonsaturable core, 24-25 
in saturable core, 24-25 
problem of, 24 
soft state, effect of, 25 

Flux division analysis, 25-48 
assumptions, 25-27 
basic equations, 28 

calculation, general, 27-32, 125- 129 

calculation, method of, 27 
circuit, 26 

computation procedure, 45 
definition of, 28 
electrical-circuit analogue, 27 
H variation in each leg, 41-43 

// in each leg, 44-45 

NT regions vs. linearity of switching 
parameters, 44-45 
parameters, 31 
summary, graphical, 37 
4>(F) variations in each leg, 30 
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Flux division, analytical expressions for: 
D - 33 

f eneral, 32, 45 

eg- dimensions effect on, 37, 45-47 
NT - 33-34, 46 

unloaded core, 32-33, 46 
8 - 1, 34- 36, 47 

Flux division in nonsaturable core, experi 
mental, 76-80 

comparison with Laddie experiment, 81 

discussion, 77, 81-82 

maximum D, 81-82 

measurement accuracy, 81 

oacil 1 ograms , 79-80 

procedure, 76-77 

results, 78-79 

Flux division in saturable core, experi- 
mental, 83-88 
discussion, 87-90 
effect of load, 88 
leg parameters, 85-86 
maximum D, 88-89 
oscil logr ants , 85 
procedure, 83 
results, 84, 88 

Flux division in saturable core, compu- 
tation of: 

effect of linearity of leg model, 90-91 
effect of load, 88, 90, 92 
effect of leg-length ratio, 93-95 
effect of t 90-91 

for 8 - 1, 90, 92 
procedure, 45, 87 

Flux measurement: 

second-order integrator, 133-134 
flux reference, 137-141 

Flux reference, 137-141 

Flux switching: 
blocking, 71-72 
closed path for, 71 
models, inelastic, see Inelastic 
switching models 


Geometry of leg: 
const ant-wi dth , 2 
effect on inel astic- switching 
parameters, 40-41 
tapered, 12 


Hyperbolic model for B{H ) , 3-4 


Inel astic- swi tching models: 
average values of, 17-24 
effect of geometry on, 40-41 
effect of <p on, 14-16 

linear B (//) - 37 
P 

linear ^ (F) , 39 

nonlinear (//) , 38-39 
P 


nonlinear 4> {F) , 40 

p 

of constant- width leg, 8 

of material, 8, 37-39 
of tapered leg, 14 . 

of tapered leg, calculation of, 14-17 
parameters of, 8, 37-41 

parameters of, experimental determination 

of, 73, 76, 151-152 
under low field, 38-41 

Integrator, second order, 133-134 


Junction problem, 60 


Laddie, 24 

flux division in, 81 

Leakage dc current in unsetting-effect 
experiment , 99- 10 1 

Leg: 

constant width, see Const ant- width leg 
dimensions of, see Dimensions, legs 
geometry of, see Geometry, leg 
lengths, 61-63 ... 

lengths, effect on flux division, 45-47, 
93-95 

parameters of, see Parameters, leg 
static 4>{F) curve of, see Static 4>(F) curve 
tapered, see Tapered leg 

Leg model limitations 
general, 7-9 

of variable- width leg, 10 

Linear material, 11 

Linear region of: 

B (H), 37 
P 

i p (F), 39 

Load in flux division: 
circuit, 25-26 
effect of, 88, 90, 92 


Machining of multipath cores, 60 
MAD: 

saturable, 25 

static 4>(F) of, 53-54 

unsetting- ef feet in, 96-101 

Magnetization switching in re-entrant 
core, 65-68 


Material parameters: 

experimental determination 

72-74 

of inelastic B (//) linear 


of, 69-70, 
region, 37-38 


of inelastic B {H) nonlinear region, 38-39 
P 

of static BiH) curve, 3-4 
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Model ; 

leg, limitations of, 7-9 
parabolic, for £, 8, 37-39 
parabolic, for p, 8, 39-40 
static B{H) , hyperbolic, 3-4 
static B(H) , linearized, 1-2 
static 0(F) 2-3, 5-6, 109-110 

vari able- width leg, limitation of, 10 

Modified parabolic model for slow 
swi tching: 
for material, 38-39 
for leg, 40, 55-56 


Nonlinear region of: 

B (H), 38-39 
b P 

<t> p (F ) , 40 

Nonsaturable core, 24-25 
dimensions of, 59 
soft state in, 25 


Osci 1 1 ograms : 

first and second order integration, 134 ■ 

flux division, nonsaturable core, 79-80 
flux division, saturable core, 85 
unsetting effect, 98 


Parabolic model for: 

E, 8, 37 

5, for low //, 38- 39 
P, 8, 14- 15, 39, 55-56 
p, for low F, 40, 55-56 

Parameters, leg, 70 

effect of geometry on, 40-41 
experimental determination of, 70-71. 

85-86, 147-152 
linear region of, 39 
nonlinear region of, 40 
test procedure, 147- 152 

Parameters, material, see Material 
parameters 

Path for leg-parameters experiments, 71 

Peak inelastic-switching coefficient: 

P . 8, 14-15 

P tp , H-15 


Re-entrant shaped core, 61, 64 

model for slow switching in, 61-69 

Reluctance, 11 

Reluctance of toroid, 115 

Reluctance of vari able- width leg: 
general, analytical, 113-115 
general, graphical, 118 
tapered leg, 115-116 
Transfluxor, 116-118 


Saturable core, 24-25 
dimensions of, 59 

Saturation flux for: 

inel asti c- swi tching model, 8 
static 0(F) model, 6 

Saturation flux density for: 
inelastic- swi tching model, 8 
static B(H) model, 3 

Soft state: 

effect on flux division, 25 
in nonsaturable core, 25 
in variable-width leg, 10 

Squares, curvilinear, 12 

Static B(//) curve, hyperbolic model for, 3-5 
coercive H, 5 
for H > H th , 4 

in saturation, 3 
material parameters, 3-4 
threshold H, 4 

Static B(H ) curve, linear model of, 1-2 

Static B(H) loop, 1 

Static 0(F) curve, experimental; 
discussion, 73 

effect of switching rate, 51-53 
procedure, 72 
results, 72-75 

Static 0(F) curve model, calculated from: 
hyperbolic B{H ) , 5-6, 109-110 
linearized B(H) , 2-3 

Switching coefficient, 8 

Switching, inelastic model, see Inelastip 
switching model 

Switching line, 61, 64, 68 

Switching path, 61, 64, 68 


Tapered leg, 7, 12 

curvilinear squares in, 12 
geometry of, l2 
reluctance of, 115-116 
P t , calculation of, 11-17 

P t , expression for, 14 
P tp , H-15 

P tp * effect of r^/r a on, 15 
P t (4 y), 22-23, 121-122 

Tests of: 

disks, blank, 53-55 
leg parameters, 147-152 
static 0(F), 149-151 

0(F), 151-152 

Threshold H, 4 
Toroid: 

properties of, 51-53 
reluctance of, 115 

Transfluxor, reluctance of, 116-118 
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Unsetting effect: 
description of, 95 
effect of fall time on, 101-102 

Unsetting effect, experimental: 

discussion, 97-102 

effect of leakage direct current in 
driver, 99- 10 1 
oscillograms, 98 
procedure, 96-97 
results, 97-99 
tested core, 96 
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